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Abstract 

We consider the Laguerre Unitary Ensemble (aka, Wishart Ensemble) of sample covariance matrices 
A = XX*, where X is an N X n matrix with iid standard complex normal entries. Under the scaling 
n = N + [\fAcN J, c > 0 and N —> oo, we show that the rescaled fluctuations of the smallest eigenvalue, 
largest eigenvalue and condition number of the matrices A are all given by the Tracy-Widom distribution 
(/? = 2). This scaling is motivated by the study of the solution of the equation Ax = b using the conjugate 
gradient algorithm, in the case that A and b are random: For such a scaling the fluctuations of the halting 
time for the algorithm are empirically seen to be universal. 
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1 Introduction 

Consider the sample covariance matrix A = XX* where X is an N x n matrix with iid entries with some 
distribution F. Our main result is a limit theorem (see Theorem |1.3[ ) for the condition number of these 
matrices wherEf ~ X Cl X c has the standard complex normal distributiorg] and 

a := n — N = [VfcN\ + o{N * 1 P), c > 0. (1.1) 

The study of this scaling is motivated by a computational problem discussed below. In the case that F ~ X c 
we refer to the matrices {A} as lying in the Laguerre Unitary Ensemble (LUE). 

°Email: deift@cims.nyu.edu, govind_menon@brown.edu, trogdon@cims.nyu.edu (corresponding author). 

1 The notation F ~ X means that the random variable F has the same distribution as the random variable X. 

-The random variable X c = X\ + 1X2 where X 1, X2 are independent Gaussians with mean zero and variance 1/2 is said to 
have the standard complex normal distribution. 
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Main results. Define the Airy kernel 

Ai(x)Ai'(y)~Ai(y)Ai'(x) 

K Ai (x,y) =-—--1 

where Ai denotes the Airy function [2D.. Then, define the Fredholm determinant 


T* 2 (s) det(7 ATai | L 2 ((s,oo))) ■ 

Here AAi|i 2 (( s ,oo)) represents the integral operator with kernel Km acting on L 2 ((s, 00 )). The function F^s) 
is the distribution function for the largest eigenvalue of a random Hermitian matrix in the edge-scaling limit 
as TV — > 00 and is known as the Tracy-Widom (ft = 2) distribution [ 251 . For a positive Hermitian matrix 
A, let A max be the largest eigenvalue of A, A m i n be the smallest and let re = A max /A m i n be the condition 
number. 

Fix c > 0. We prove: 

Theorem 1.1 (Smallest eigenvalue A m i n ). Assume A = XX* where X is an N x (TV + a) matrix of iid 
standard complex Gaussians and a = \_V 4cTVJ . Then for all t eR 


lim P 

N—too 


C ^min 
ca -2 / 3 2 2 / 3 



F-2(t). 


Theorem 1.2 (Largest eigenvalue A max ). Assume A = XX* where X is an N x (N + a) matrix of iid 
standard complex Gaussians, a = [ \/ 4c and v = 4iV + 2a + 2. Then for all t 


lim P 

N—>oo 


( ^max V 

U 1 / 3 2 2 /3 



F 2 {t). 


Theorem 1.3 (Condition number re). Assume A = XX* where X is an N x (TV + a) matrix of iid standard 
complex Gaussians, a = \pjAcN\, and v = 4TV + 2a + 2. Then for all t 6l 


lim P 

N—> 00 



z/a _ 2 / 3 2 2 / 3 



lim P 

N—>oo 


4 N 


4c —4 / 3 A^ 2 / 3 



F 2 (t). 


History of the problem. The study of the eigenvalues, and in particular, the condition number, of 
random positive definite matrices has a rich history in mathematics and statistics going back at least to 
the seminal paper of Goldstine and von Neumann m- The exact distributions of the largest and smallest 
eigenvalues of sample covariance matrices, with iid columns, were computed in |2.'ll and El, respectively, in 
terms of infinite series and hypergeometric functions for any finite N and n. When F is either a standard 
real or standard complex Gaussian distribution and a = 0, Edelman [7] determined the scaling limit of the 
smallest and largest eigenvalues and the condition number as N —> 00 , 

p (jP Sf ) = FB « : = e_4, ‘- < L2 > 

It also can be shown that the condition number distribution is heavy-tailed for finite N because the density 
of A m i n does not vanish near zero — A m ; n is exponentially distributed with parameter N/2 [?]. 

As noted by Johnstone]^] [15], a by-product of Johansson’s work on last-passage percolation [14] is that the 
fluctuations of the largest eigenvalue of LUE with a fixed are given in terms of the Tracy-Widom distribution 
7*2(t) as N —> 00 . When a = [_c7VJ, c > 0, it can be shown that the smallest eigenvalue also has Tracy- 
Widom fluctuations in the TV —> 00 limit [T| and that the condition number has fluctuations given by the 
convolution of two independent, but scaled, Tracy-Widom distributions. This result does not appear to be 
explicitly stated in the literature but it follows from the asymptotic independence of the extreme eigenvalues 
[2]. See jT3] for the case of n> TV 3 . 

3 Johnstone m also proved a limit theorem when F is a standard real Gaussian. 
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In terms of interpolation between the limiting condition number distribution F^(t) at a = 0 and the con¬ 
volution of two Tracy-Widom distributions when a = |_ciVJ, we see that the scaling (1.1) is sufficiently strong 
to force A m in away from zero and give pure Tracy-Widom statistics for the condition number. Assuming 
independence of the smallest and largest eigenvalues, from Theorems o and| 1 . 2 | we have 


An 


V + i/ 1 / 3 2 2 / 3 ^Qy 


GUE 


Amin c — Ca _2 / 3 2 2 / 3 ^Q 1 l 


GUE 


Here £qu E and £que are 'FI random variables with P(£gue — b) = -^(t). Then using a Neumann series we 
have a formal expansion 


t(2) 


( 1 ) 


An 


An 


V 

c 


(»- 2 «2n<?i E + . 


(1.3) 


If a = |_ciVJ, the first two terms in this expansion are of the same order and dominate the expansion. Thus, 
it is clear why the convolution is involved. Also, for a = 0 it is clear that such an expansion cannot be 


justified. In the case of (1.1), a. <C v and just the first term is dominant. Lemma 4.6 rigorously justifies the 


formal expansion in this case. It does not appear that the scaling (1.1) has been treated previously in the 
literature. 

Remark 1.1. Our calculations in this paper are for the case a = [(4c./V) 7 J, 7 = 1/2. Note that the first 
term in the right-hand side of (1.3) is still dominant for 0 < 7 < 1. For this reason we anticipate a similar 
Tracy-Widom limit theorem for the condition number for 0 < 7 < 1. In addition, we anticipate that the 
conclusions of the numerical studies that we discuss below when 7 = 1/2 will be unchanged for 0 < 7 < 1 . 

Our method of proof makes use of the asymptotics of Laguerre polynomials x)} n >o which are 

orthogonal with respect to the weight x a e~ x dx on [0, 00 ). These asymptotics are derived via the Deift-Zhou 
method of nonlinear steepest descent as applied to orthogonal polynomials [4] (see also [3] for an introduc¬ 
tion). For fixed a , this problem was addressed by Vanlessen [27| for the generalized weight x a e~^^dx (see 
also [2D for Q(x) = x). From the classical work of Szego [24] and [27], the asymptotic expansion of L^\x) 
as N —>■ 00 near x = 0 is given in terms of Bessel functions. Forrester [9] noted that this implies the statistics 
of the smallest eigenvalue are given in terms of a determinant involving the so-called Bessel kernel. 

As is seen in Theorem under the scaling ( |1.1[ ) , the asymptotics of Laguerre polynomials is given in 
terms of the Airy function, giving rise to the Airy kernel and producing Tracy-Widom statistics. This was 
noted first by Forrester [5]. This transition from Bessel to Airy can also be seen by considering the weight 
x a e~ x ~ t / x dx for varying t [28]. The difference here is that this transition is induced via the parameter a 
that is naturally present in the Laguerre polynomials. 


A computational problem. Our motivation for considering the scaling 0 comes from a computational 
problem. In numerical analysis, the condition number of a positive-definite N x N matrix A is arguably the 
most important scalar quantity associated to the matrix. Specifically, it controls the loss in precision that is 
expected when solving the system Ax = b. The condition number can also be tied directly to the difficultly 
encountered in solving the system by iterative methods. This is evident, in particular, in the conjugate 
gradient algorithm used to solve Ax = bm i- The conjugate gradient algorithm is stated as follows (see E3 
for an overview): Given an initial guess xq (we use Xq = 6 ), compute ro = b — Axo and set po = r o- For 
k = l,...,N, 


1 . Compute the residual 77 = 77_1 — a k -\Ap k _\ where 


«fc- 1 = 


(r fe i, r fc i) 


e 2 


(Pk-i,Apk-i)p 


2. Compute p k = ?7 + b k -iPk-i where b k -1 = 


3. Compute x k = x k _i + a fe _ip fc _i. 


{r k ,r k )p 
(r k ~i,r k -i)p ' 


4 We use the notation I 


11^2 


= (v,y)e 2 = E,- l%j 2 for y = ( 2 / 1 , 272 ,... ,vn) e c N . 
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If A is strictly positive definite Xk -+ x = A l b as k —> oo. Geometrically, the iterates Xk are the best 
approximations of x over larger and larger affine Krylov subspaces ICk, 

|| Ax k — b\\ A = mm xGK : k \\Ax - b\\ A , IC k = £o + span{r 0 , Ar 0 , .. ., A fc_1 r 0 }, ||x||^ = {x,A~ 1 x) i 2 , 

as k f -IV. In exact arithmetic, the method takes at most N steps: In calculations with finite-precision 
arithmetic the number of steps can be much larger than N. The quantity one monitors over the course of 
the conjugate gradient algorithm is the norm ||T~fc ||^ 2 . By m , 


\\rk\\e 2 <2(1 + -J= 


-2k 


Folk 2 , 


(1.4) 


from which we see that the larger k is, the slower the convergence. It is remarkable that the bound (1.4) 


does not depend explicitly on n, only implicitly through k. We note that (1.4) was derived in m under 
the assumption of exact arithmetic, but (1.4) is also useful in calculations with finite precision, provided the 


effect of rounding errors is anticipated to be small. 

We are interested in the statistical behavior of the conjugate gradient algorithm (l)-(3) when A and b are 
chosen randomly. Let b have iid entries distributed according to a distribution F which may differ from the 
matrix-entry distribution F. We use the pair E = (F, F) to refer to the ensemble, encoding the distribution 
of the entries of both A and b in Ax = b. Let e > 0, E, N > 0 and n > 0 be given. For a pair (A, b ) we define 
the halting time T e ,E,N,n = T^e,n,u{A, b ) to be the smallest integer such that ||r-fe ||^2 = ||rfc(A, 6)||^2 < e. The 
residuals in the conjugate gradient method decrease monotonically, thus ||rfc(A, 6 ) ||^2 < e for k > T ei E,N,n- 

In [6], the authors performed a numerical study of T e ,E,N,n as a random variable. For n = N + [2%/iVJ, 
Monte Carlo simulations for different ensembles E were used to show that the fluctuations 


T e ,E,N,n 


T e ,E,N,n ~ lE[T e ,g,JV,w] 
\/Var[ T e ^E,N,n\ 


of T e E,N,n had a universal limiting form 5 1 The motivational goal of this paper is to further our understanding 
of the universality that is observed in jo] under the scaling (1.1). To this end, we focus on the condition 
number n = = kn,u{A) of the matrix A and use the the rigorous results above, numerical simulations 

of T tj E,N,n 33 and Kjv.n, and the estimate (|1.4|), to infer properties of the halting time distribution. 

Numerical results presented in Appendices |A.l| |A.3| illustrate the sharply different behavior of the condi¬ 
tion number and halting time distributions in the cases a = 0, a = [2v / )VJ and a = N. Universality appears 
to depend critically on the choice a = \_2y/N\: If a = 0 or a = N, the histogram for T e E,N,n does not 
appear to have a universal form. The numerical experiments also reveal an interplay between the tightness 
of the condition number distribution and tightness of the halting time distribution that is consistent with 


the following upper bound on T e E,N,n , which is obtained by taking a logarithm of (1.4): 


1 log 11r 0 ||.£2 - loge + log2 

-L e.E,N,n S 


l0 § {} + 7s) 

In order to use this bound to estimate the expectation of T^ E N n , j = 1,2,... one must estimate 


E 


We make some observations: 


lo s J i 1 + 7s) 


lo S J i 1 + 7s) 


dP(/c < s). 


(1.5) 


(i) Since log- 7 ^1 + -^j = 0(s H 2 ) as s —» 00 , (1.5) does not provide an a priori bound on the higher 
moments of T t ^E,N,m N and n fixed, if the distribution P(kjv j71 < s) is heavy-tailed. 


(ii) The right hand side of (1.51 may diverge as N —> 00 even when all moments are finite for finite N. 


3 By universal, we mean that the histogram for r e ^E,N,n is independent of E for N sufficiently large and e sufficiently small. 
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(iii) When the system is well-conditioned, for example if kn, u — C+N ~ a £, a > 0, C > 1, for a fixed, positive 
random variable £ with exponential tails, all moments of T €t E,N,n are bounded by an TV-independent 
constant. 



This is consistent with the numerical results, since the empirical mean of the halting time is (much) larger 
than TV. Since the maximum number of steps in the conjugate gradient method is TV in exact arithmetic, 
this also shows that round-off errors have degraded the accuracy o f the computation. 

The well-conditioned case, n = 2TV, is considered in Appendix A.2 We assume F is uniform on [—1,1] 
and F ~ X c or F is Bernoulli. In contrast with the case above, the condition number Kj\r n satisfies (iii) in the 
limit TV — > oo with a = 2/3 implying, in particular, by (1.5 1, that supjy E[T £i B,JV,n], sup^ Var[T e ^,JV,n] < oo. 
The numerical results in Figures 12(b) and 12(d) indicate that the random variable T tt E,N,n remains discrete 


in the large TV limit. Indeed, this is necessarily true for any sequence of integer-valued random variables 
(A/v)at>o with supjy VarJTX/v] < oo: Clearly, in such a case the fluctuations 

X N — E[A/y] 

VVar[A/v] ’ 

cannot have a limiting distribution with a density. The only universal limit possible here is a statistically 
trivial point mass. 

Finally, we turn to the critical scaling, n = |_2-\/iVJ in Appendix 


A.3 


It follows from Theorem 


1.3 


that 


E[«jv,n] = 0(N), Var [n N , n ] = 0(N 4 / 3 ), TV -> oo. 


As in (ii), the estimates on the right hand side of (1.5) diverge as TV —»• oo. This divergence is consistent 
with an empirically observed divergence in the mean and variance of T t ^E,N,n- Indeed, in contrast with the 
well-conditioned case, such a divergence is necessary to obtain a non-trivial limiting distribution for r e E,N,n- 
Further, in contrast with the ill-conditioned case, the empirical mean of the halting time, while large, is 
(much) smaller than TV and rounding errors do not appear to play a dominant role. 

The paper is organized as follows. In Section [2] we review the global eigenvalue density for LUE and 
discuss its connection to Laguerre polynomials and hence to a Riemann-Hilbert problem. In Section [3j we 
use classical Riemann-Hilbert analysis to rigorously determine the asymptotics of the Laguerre polynomials 
that appear in the global eigenvalue density. In Section [4] we use these asymptotics to prove limit theorems 
for the distribution of the largest and smallest eigenvalues, along with the condition number. This final 
section contains the proofs of the main results. As the behavior of the conjugate gradient algorithm is 
universal with respect to the choice of E above, it is sufficient to consider one particular ensemble. For this 
reason we have decided to study the analytically tractable case E = LUE. We include a table of notation to 
guide the reader. 


Acknowledgments. This work was supported in part by grants NSF-DMS-1300965 (PD), NSF-DMS- 
1411278 (GM) and NSF-DMS-1303018 (TT). The authors thank Ivan Corwin, Anne Greenbaum, Ken 
McLaughlin and Brian Rider for helpful conversations. 


Nomenclature 

Ai The Airy function, page 1 
a a = N — n, page 1 

w w(z) = (2a + 2 )z — a log^ z + (a + l)7ri, page 22 
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t N l N = 27V(21og2 + 1), page 12 

w w(x) = (2a + 2)x — a log a:, page 12 

H))- 1 . H^ Hankel functions, page 18 

I Q , K« Modified Bessel functions, page 18 

log < _(-) The principal branch of the logarithm, page 10 

(')Z- The principal branch of the root, page 10 

/Cn The correlation kernel, page 7 

AT The solution of the RH problem with a skew jump, page 15 

v v = 41V + 2a + 2, page 2 

V The left factor for the correlation kernel near z = 1, page 31 

W The right factor for the correlation kernel near z — 1, page 31 
</><_ The indefinite integral of the EM cut along (— oo, 1], page 10 

ct !>_► The indefinite integral of the EM cut along [0, oo), page 10 

Vv = {zfj 2 + {z- l) 1 / 2 , page 12 

i>^(z) = (z)^ 2 + (z - l)^ 2 , page 11 

log^(-) The branch of the logarithm that is cut on [0, oo) with a real limit from above, page 10 
(■)\ The branch of the root that is cut on [0, oo) with a positive limit from above, page 10 
T e ,E,N,n The halting time fluctuations, page 4 

A The solution of a Rienrann-Hilbert problem near z = 0, page 24 

c The constant in the critical scaling a = |_%/4ciVJ, page 1 

D The solution of a diagonal RH problem, page 16 

Coo D( oo), depends on IV,a, page 16 

/«_ The change of variables near z = 1, page 20 

/_► The change of variables near z = 0, page 22 

g The log transform of the EM, page 10 

M M = N + |(a + 1), page 20 

Ma; The analytic prefactor for S<-(z), page 20 

N The dimension, N x N, of the LUE, page 1 

PAi The classical Airy parametrix, page 17 

Peessei The classical Bessel parametrix, page 18 

S'oo The solution of the truncated RH problem on [0,1], page 15 

S The local parametrix near 2=1, page 20 
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S'-). The local parametrix near z = 0, page 22 
T e ,E,N,n The halting time, page 4 

V The left factor for the correlation kernel near 2 = 0, page 29 
W The right factor for the correlation kernel near 2 = 0, page 29 
w v w v (x) = x a e~ l ' x , page 8 

EM The Marchenko-Pastur equilibrium measure, page 10 
LUE Laguerre Unitary Ensemble, page 1 

2 The Laguerre Polynomials and the Laguerre Unitary Ensemble 

2.1 The Laguerre kernel 

We first recall the definition of the Laguerre kernel. These ideas are well-known and may be found in 0 
Section 2]. Let A = XX* where X is an TV x (N + a) matrix of iid standard complex Gaussian random 
variables. Then the eigenvalues 0 < A m i n = Ai<A 2 <---<Aat = A max of A have the joint probability 
density 


Pn( Ai, ..., A N ) 


1 N 

n i^»i 2 - 

°AT j —1 l<j<k<N 


The statistics of eigenvalues are more conveniently expressed as determinants involving Laguerre polynomials. 
Recall that the Laguerre polynomials, {L^ a \x)}^. 0 , are a family of orthogonal polynomials on [0, oo), 
orthogonal with respect to the weight e~ x x a . We normalize them as follows [20] 


L^\x) = kjxi + 0{x^ 1 ), kj = 


(~iy 


J ji 


/»oo 

/ L ( i a \x)L < j a \x)e~ x x a dx = Si 
Jo 


T(j + a + 1) 


( 2 . 1 ) 


We then define the associated wavefunctions, orthogonal with respect to Lebesgue measure on [0, oo), 


4>j{x) := 


J ■ 


1/2 


r(j + a + 1) 

il)j(x)ipi{x) dx = Sij, 


e~ x/2 x a/2 L { y\x), 


and the correlation kernel 


N-l 

^n{x,v) = y; il>j{x)ipj(y), 0 < x,y < oo. 
j =o 

The kernel /Cjv defines a positive, trace-class operator on L 2 ([a,b\). Since has finite rank, it is clearly 
trace class. To see that ICn is positive, consider / G C°°((s,t)) with compact support and note that 


n t ft ntN— l 

JC N (x,y)f(x)f*(y)dxdy = / / '^2 l il)j{x)ip j (y)f(x)f*(y)dxdy 

J S J S j—Q 


I yyx)f{x)dx S j (^1 ip j (x)f*(x)dx S j = 


ipj(x)f(x)dx 
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It is by now classical that the statistics of the eigenvalues Ai ... < Xn may be expressed in terms of 
Fredholm determinants of the kernel /Cjv0[3] (which are well-defined since KLn is trace-class). In particular, 
the statistics of the extreme eigenvalues are recovered from the determinantal formula 


P (no eigenvalues in [a, 6]) = det(I — /Cjv| L 2 ([ a ,b]))- 
By the Christoffel- Darboux formula [M] , we may also write 

N\ {T(N + a + 1)T(N + a)\ 1/2 (i/; N - 1 (x)ip N (y) - tp N (x)ip N ^ 1 (y) 


( 2 . 2 ) 


K, N {x,y) = 


T{N + a) 
N\ 


N\(N — 1)! 


x-y 


r -(x+y)/2 a/2 a/2 L N ^ (v) L N-l( X ) ^ (af^AT-l (v) 

T(N + a ) V x-y 


Thus, questions about the asymptotic behavior of K,N{x,y) as N 
asymptotics of ant 
regime for a, see (1.1). 


oo reduce to the study of the large N 
asymptotics of and What is new in this paper is the study of these asymptotics in the scaling 


2.2 The Riemann-Hilbert approach to Laguerre polynomials 


To compute the asymptotics of the Laguerre polynomials we use their representation in terms of the solution 
of a Riemann-Hilbert problem and follow m for the general theory. We also refer to [3Tj for some explicit 
calculations. We note that we use Riemann-Hilbert theory as opposed to using the integral representation 
for Laguerre polynomials to determine the appropriate asymptotics. This is because in the scaling region of 
interest the Riemann-Hilbert method gives a direct and algorithmic approach to the difficulties arising from 
turning point considerations. Define the rescaled polynomials (cf. CHI ) 

L^\x)=Lf\vx), v = 41V + 2a + 2, (2.3) 

n j( x ) = (— 1 Y j\v~^ L^\x) = x n + 0(a; n_1 ). (2.4) 

Here { 7 ^ ( 2 ;)} are the rnonic orthogonal polynomials for the weight 

w v (x) = /e" ra . (2.5) 

This scaling is chosen so that the asymptotic density of the zeros of Lff (x) as N —> 00 is supported on the 
interval [0,1] (see S21] for a fixed). Following [8] we define 


Y(z) = 

where Cr denotes the Cauchy integral operator 


tt n (z) C r+ [ir N w„} (z) 

c N ir N -i{z) c n C r +[tt N -iw„](z) 


, zeC\l 


and, by (2.1), 


c" 1 - - 
N ~ 2 TTi 




( 2 . 6 ) 


In the remainder of the manuscript we use the notatiorj^] Y ± (z) = Y±(z ), z £ T to denote the boundary 
values of an analytic function Y(z) from the left (—) or the right (+) side of as one moves along an oriented 
contour T. See Figure [I] for schematic of a contour F. We also use the notation 

Cpf(z) :=(C r /(z)) ± . 


By general theory, the matrix Y (z) 
Riemann-Hilbert problem: 


defined in equation (2.6), is the unique solution to the following 


We allow the ± to be in either the sub- or super-script for notational convenience. 















Riemann Hilbert Problem 2.1. The function Y{z) satisfies the following properties: 

1. Y(z) is analytic on C \ [0, oo) with limits Y ± (x ) as z —> x £ (0, oo) from above or below; 


2. Y+(x) = Y~(x) 


1 w„{x) 
0 1 


for x £ (0, oo); and 


3. Y(z) = (/ + 0(z *)) z N<73 as z —> oo. Here a 3 = diag(l, —1). 

We focus our effort on the analysis of Y(z), since this yields the asymptotics of K-n via the identity 
k- , , w*(x)wkv)[ 0 1 ] Y+\x)Y + (y) [ 1 0 ] T ^ 

IC N {x : y) = -—----- k , 0 < x,y < 00 . (2.7) 


27ri 


x-y 


Indeed, the definitions of the limits Y ± , the weight w u , and the rescaled Laguerre polynomials (2.3), yield 
the identity 


)C N (x, y) = wf (x)wS {y) 


N\v a L^iy)^ n-i(x) - L ( * ] (x)L^li(a;) 


F(N + a) 


= w3 (x)w3 (y)zzTW - r(v~ 2N+1 )N\ 


x-y 

Lm\v)Lm- iW - (xjL^l^x) 


r (N + a) 

1 1 

w3 ( x)w3 (y) n N (y)n N _ 1 (x) - TT N (x)Tr N -i(y) 


x-y 


1 

c N 


(27Ti) 


x-y 


w3 ( x)w3 (y) Y+{y)Y+{x) - Y+(x)Y+(y) 


27ri 


x-y 


Since det(P( 2 :)) = 1, this identity is equivalent to equation (2.7). 



Figure 1: An example of an oriented contour T. 

The goal of the reminder of this section is to solve for Y(z) asymptotically as N — >• 00 . This requires a 
series of explicit transformations or deformations: 

• Y 1 —> T so that T(z) ~ / as z —> 00 , 

• T > S so that the jump matrices for S tend uniformly to the identity matrix on closed subsets of 
C \ [0,1], and 

• S 1 —> E where the jump matrices of £ tend to the identity matrix in L 2 D L°°. 

This procedure is now standard and a general reference is |2j. 
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2.3 The first deformation: normalization at infinity 

Before we proceed we introduce additional notation to fix branch cuts. For 7 £ R. let 

Z *( Z )l=z>, 

have its branch cut on (— 00 ,0] such that (z)7_ > 0 for z > 0, i.e., the principal branch. By contrast, let 

* -> {z )\. 

have its branch cut on [ 0 , 00 ) in order that (z)Z> + > 0 for 2 > 0 , when the real axis is oriented left-to-right, 
i.e. the limit from above is positive. Note that for Im z > 0, {zfU = If 2 < 0, then (z)\ = |z|e 1 ’ 1 ' 7 . 

We similarly define branches of the logarithm. The map 

z !->• log^(z) = log(z), 


denotes the principal branch of the logarithm, whereas 


2 ^ log_>(z), 


has its branch cut on [ 0 , 00 ) so that log_^ + (z) > 0 for 2 > 1 . 


It is then clear that 


l°g_>.(z) 


log<_(;z), if Imz > 0 , 

log.,_(;z) + 27ri, if Imz < 0. 


In this notation the left arrow is used for emphasis and signifies the principal branch. It is sometime 
omitted when there is no confusion. Also \f■ is always used to denote the non-negative square root of a 
non-negative number. 

As in [5j, we remove the polynomial behavior of Y ( z ) at infinity using the log transform of the equilibrium 
measure as follows. In our case a/N -> 0 as A -> 00 and the equilibrium measure (EM) is the so-called 
Marchenko-Pastur distribution m 


Ms) = ~ ( s )ds, 

where xa denotes the characteristic function of the set A. We define the log transform 

g(z) ■■= [ log_,(z- s)dii(s). 

J 0 


( 2 . 8 ) 


(2.9) 


We also introduce the following functions to simplify the analysis of the asymptotic behavior of g as z —> 00 
and its jumps across the interval ( 0 , 00 ): 


f z ( s - 1 \ 1/2 

j ds, z€C\[0,oo), 

(j> i -(z):=2j ^ ds, z£C\(-oo, 1]. 


( 2 . 10 ) 

( 2 . 11 ) 


Lemma 2.1. (i) For x > 0, the upper and lower limits of g are given by 


g + {x)= log \x - s|d/Lt(s) +7ri / d //(s), 

Jo Jx 

g~(x) = / log \x — s|d/i(s) + 27ri — 7ri / d/i(s), 

Jo Jx 

9 + {x) - 9~(x) = / d n(s). 

Jo 
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(ii) e Ng( ~ z l = z N { 1 + C>( 2 _1 )) as 2 —» oo, uniformly in C \ [0, oo). 

(in) For z > 0 

g ± (z) = -<j)%(z) +2z - (2 log2 + 1) + 7 ri, (2.12) 

4>%{z) = ™ + <A<-(z), (2-13) 

<t>~f{z) = —i 7 r + ((<-{z). (2.14) 

Proof. Parts (i) and (ii) follow directly from the definition of g(z) and log_>.( 2 ). To establish (iii), we first 
show that 


g(z) + <t>^(z) 


extends to an entire function of z. First, it is clear that the function is bounded in the finite plane, and 
analytic for Re z < 0. For Re z > 0 we check the boundary values. For 0 < Re z < 1 we have 


9 + {z) + <t>%{z) =g (z) — 2ni f dg(s) + </>+ (z) = g (z) + (j>^(z). 

Jo 


(2.15) 


Now, if 2 > 1 we use (2.13) and (2.14) and (2.15) holds. This shows that g(z) + 4>^{z) is an entire function 
because 0 and 1 would have to be bounded, isolated singularities. To establish (2.12) we turn to an explicit 
integration of <^( 2 ). Consider for 0 < z < 1 


nz 1 1 _ 2 py/z 

</>+( 2 ) = 2i J y-ds s = 4i J \/l — t 2 dt = 2i(y / z{ 1 — 2 ) + arcsin y/z). 

Then, we must consider the analytic continuation of this function so that it is analytic in C \ [0, 00 ). We find 
that if 6 = arcsin y/z then 

0 = -ilog < _(iy / i ± \/l — z 2 ). 


If we take the (—) sign, we actually need n — 6 to consider the correct branch. From this, it follows that 
2 i(i/ 2(1 — z ) + arcsin y/z) = 2 (2(2 — l))j( 2 + + 27ri — 2 log^(i ( z ) 1 ^ 2 + i (2 — 1 )^ 2 + ). 

Define 

iM*) = (z)T + (z 1 )]i 2 • 

First, it can be shown that if if^-(z) = c £ R. then 2 > 1. But for 2 > 1, it ip~y(z) > 0 so that must map 
into either the upper- or lower-half plane. We find that Imip^(z) > 0 for 2 £ C \ [ 0 ,00). From this we find 
the analytic continuation to the whole complex plane 

< j >- t ( z ) = 2 (2(2 — 1)) 1,/2 - 2 ^^ if - tiz ) + in, z G C \ [ 0 , 00), 


where log^ an be replaced with log^ because lm , ij;_ ¥ (z) > 0. Also, for this reason, log^ i/j^.(z) is analytic 
anywhere if)- ¥ (z') is. We also have 


(f^(z) = 2z — log 2 — (2log2 + 1) + 7 ri + o(l), g(z) = log 2 + o(l), as 2 —> 00 . 


This proves (2.12). 

We now use this proof to obtain a closed-form expression for <f>^_{z). First, </>_>. ( 2 ) — 7 ri 
Im 2 > 0 and then for 0 < 2 < 1 


<£<-( 2 ) = 2 ( 2(2 - 1))+ /2 - 21og^(i(2)^ 2 + + i (2 — 1)(^ 2 + ) + 7ri 
= 2 ( 2(2 - 1))+ /2 - 2 log^_(i ( z)]!_ 2 + + i (2 - 1)^ 2 +) + 7ri 
= 2 ( 2 ( 2 - l))+ / 2 - 21 og^ ^+(z), 

VU z) = (z ) 1 / 2 + (z-l ) 1 / 2 . 


(2.16) 

□ 

for 
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(2.17) 

(2.18) 













Therefore, the analytic continuation is given by 

= 2 (z(z - 1)) 1/2 — 2 log^ ip<_(z), z £ C\ (— 00 ,1]. 

and > 0 for z £ C \ (— 00 ,1]. 

Define = 2TV(21og2 + 1) and the matrix function 

T(z) := e^ eNa3 Y(z)e-( N9( ' z ' >+ 5 eN ^ a3 , 

which solves the following problem. 

Riemann Hilbert Problem 2.2. The function T{z) satisfies the following properties: 

1. T(z) is analytic in C \ [0, 00 ) and has limits T ± (x) as z —> x £ [0, 00 ) from above or below. 

2. T + (x) = T~(x)Jt(x) for x £ [0, 00 ) where 

p -w(x) 


Jt(x) = < 


0 


e 

p 2N<fi~ t (x) 


p—2,N(j>~fi(x)—w(x) 

0 1 

w(x) = ( 2 a + 2)x — a log x. 


if x £ (0,1), 

, if x£ [l,oo). 


3. T(z ) = I + 0(z *) as z 00 . 


(2.19) 


2.4 The second deformation: lensing 

We factor the matrix J T on (0,1) as follows: 

J T (x) = 


1 0 ' 


0 e -M=o) ■ 


1 

0 ' 

e 2N <t>Z,(x)+w(x) 1 


—e^C*) 0 


(x)+R>(ai) 

1 


This allows a lensing of the problem. Let T| and Tj. be contours as in Figure [2] and set 


( 2 . 20 ) 


,5(2) := < 


T(z), 

T{z) 

T(z) 


1 0 

_ e 2N(f>->(z)+w(z ) ^ 

1 0 

e 2N4>^(z)+w(z) ^ 


if 2 is outside the region enclosed by T-f and T^., 
if 2 is inside the region enclosed by [ 0 , 1 ] and Tf, 

if 2 is inside the region enclosed by [ 0 , 1 ] and Tp 


We obtain the following Riemann-Hilbert problem for S. 

Riemann—Hilbert Problem 2.3. The function S(z) satisfies the following properties: 

1. S(z) is analytic on C\T, 

2. S + (s) = S~ (s)Js(s) where Js is defined by 


Js(s) = < 


1 0 

g 2iV0_ > (s)+u)(s) 1 

0 

-e*M 0 

l e -2N4>fi(s)-w(s) 

0 1 


if S £ T t UTp 
if s £ (0,1), 
if s £ (l,oo). 
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Figure 2: The jump contours T for S. The region bounded by F^ and F| is called the “lens”. 


3. S(z ) = 0(1) as z —> 0 from outside the lens, 


l S(z) 


1 0 

± e 2N<t>^(z)+w(z) i 


= 0(1) as z —» 0 inside the lens. The (+) sign is taken for z in the 
region enclosed by [ 0 , 1 ] and T| and the (—) sign is taken in the region enclosed by [ 0 , 1 ] and T^ ; 

5. S(z ) = I + 0(z _1 ) as z —> oo. 

This deformation is valid for any choice of contours T^ and T^ arranged as in Figure [2] 


Lemma 2.2. With a as in (1.1), there exist positive constants cs and D$ such that 

I Us - ^||L“(r T \(B(o,<5)us(i,<5)))) II Js — J|U°°(r. l .\(B(o, 5 )u.B(i, 5 ))) < Dse~ Ncs , 

II-7s - jr |U 1 nL'»((i+5,cx))) < Dse~ Ncs . 

1/2 

Proof. We now obtain an upper bound on the real part of ()>_,. For Im z > 0 the function (z)' agrees with 
the principal branch. For z = x + iy, 


Im 


z-1 


y 


x 2 + y 


,2 ' 


It is also clear that for argz £ (0, 7 r), Im (z)]^ > \\z\ 1 / 2 Imz. From this, the simple estimate follows for 
Im 2 > 0 


Im 


x - 1 


1/2 


1 

> 


Im z 


2|2-l|l/2|.|3/2- 

Assume 0 < Rez < 1 , Iinz > 0 and consider the real part of the integral 


s — 1 


1/2 


ds = 


fRez / s _ | \ l/ 2 


ds ■ 


>,+ 


f (^i 

JRez \ S 


1/2 


ds. 


The first term is purely imaginary and the only contribution to the real part is from the second term. Then 


Re<^_>(;?) < — / 
Jo 


| Im z\ 


\y\ 


0 |(Re2-l) 2 +2/ 2 | 1 /4|(Re^)2 + y2|3/4 


d y <- 


r*| Im z\ 


_ \y\_ 

'o (1 + y 2 ) 


2 d y- 


( 2 . 21 ) 


This quantity on the right-hand side is bounded uniformly above by a negative constant provided Imz > 
S > 0. A similar argument follows for Imz < 0 resulting in the same estimate. 

From this discussion, it follows that 


II Js ~ -^l|z/ oo (r t \(B(0,5)UB(l J (5)))) ||- -f||l,“( r4 .\(B(0,5)015(1.,<5))) < D S e 


—Ncg 
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for some positive constants Dg,cg. For 2 >lwe also look to estimate the real part of </>_> (z). It follows that 

1/2 


Re </>+(;?) = Re <j>^L(z) = J 


r 

5—1 

h 

5 


rls, 


which is necessarily a monotonic, strictly increasing function giving the estimate 


11-7s ~ -7||i 1 nL°°((i+5,oo)) < D' s e N ° s , 

for some new positive constants D' S: c' s . These estimates hold even in the case a —> oo because e -( 2a + 2 ) z + al °s z < 
1 for 2 > 1. The lemma follows from a redefinition of eg, Dg. □. 

We require the following lemma in the sequel. 

Lemma 2.3. 


<^-►0)^0, z€C\[0, oo), 

0+(O)=O, (2.22) 

<f%.{z) ^ 0, z > 0. 


Proof. We first claim that 


signlm 


z-1 

V 2 


1/2 


= sign Im z. 


Assume that there exists two points a, 6 £ C + so that 


Im 


a — 1 


1/2 


> 0 and Im 


6-1 


1/2 


< 0 . 


(2.23) 


/ \ !/ 2 

It follows that there exists a point z* on the line that connects a and 6 so that Im ( 1 =0. Then 

(z* — 1 )/z* > 0 contradicting that z* is in the open upper-half plane. Similar considerations follow in 
C“ and the sign of the imaginary part is constant in each open half plane. Then (2.23) follows from the 

I Icy 

Cauchy-Riemann equations as (—=A) is real with a positive derivative for z < 0. 

1 /2 

Similarly, it can be shown that Re (a=T) > 0 for z £ [0,1]. We write 


/(«) = 2 Re 




1/2 


, g(s) = 21m 


s — 1 
s 


1/2 


nr\,ez f'txez 

<?M~) = / /+(s)ds + i / ff+(s)ds 

Jo Jo 

/*Imz nlraz 

+ i / /(Re 2 + is)ds — / <?(Re z + is)ds 

Jo Jo 


Then for Re 2 > 1, z G C+, 


For Rez < 1, z £ C+, 


/*Im z 

Im </>_). (z) = 7r + / /(Rez + is)ds > tt. 
Jo 


/*R ez />Imz 

Re </>_>. (z) = / / + (s)ds— / (/(Rez + is)ds < 0 

Jo Jo 
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because g is positive in the upper-half plane, f+(s) = 0 for 0 < s < 1 and f+(s) > 0 for s < 0. Also, for 
* < 0, 0-+ (z) = / 0 R °" f(s)ds < 0. Hence cannot vanish in C+ \ [0,oo) and 4>%(z) cannot vanish for 

z > 1. 

For ^ € C“, (j>-*{z) = </>_> (z*)* because the two functions are equal for z < 0. Thus, it remains to show 
that 4>%(z) vanishes on [0,1] only at z = 0. From 



Im^>_>(,£) is a strictly monotone function on ( 0 , 1 ) and the lemma follows. □ 

This leads us to consider the solution of the Riemann-Hilbert problem obtained by removing the jumps 
on F|, Tj, and [l,oo). 

Riemann—Hilbert Problem 2.4. We seek the function Soo(z) with the following properties: 

1. Soo(z) is analytic on C \ [0,1], 

. r o e -w{x) 

2- s+(x) = Soo(x) Q for x€ (0,1), 

3. S^(x) diag |x| _ 5 a ^ £ L 2 ((0,1)), and 

4- Soo(z) — > I as z —> oo. 

Remark 2.1. A standard argument shows that Soo with these properties is unique if it exists. 

We now show that exists by an explicit construction. We first start with the determination of a 
matrix function A f which satisfies the following conditions 

J\f + (x) = J\f~ (x) ^ 0 , x e (0,1), Af(oo)=I. 

Consider 

U := ~^= I - I , so that ** } = U diagli, — i)C/ _1 . 

v/2 L 1 _1 J L ~ 1 0 

Direct calculation shows that A4 = U~ l J\fU satisfies the conditions 

M + (x) = M~(x) q _0. , a; G (0,1) M(oo) = I. 
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We write the jump condition for Soo{z) as 


5+(a;) = S^ix) 


__ e w(x) 


e -w(x) 

0 


Af + 1 (x)S+(x) = 


0 -1 

1 0 


-A/L (tfS^ix) 


0 Q-™^) 

-e^) 0 


(2.24) 


We make the ansatz D(z) = diag(d(z), l/d(z)) where D(z) := M 1 ( 2 )S' 00 ( 2 ;). Then ( 2.24| ) turns out to be 
equivalent to the condition d + (z)d-(z) = e w ( z \ Formally, by taking the logarithm and letting h(z) = log d(z) 

h + (z) + h~(z) = ( 2 a + 2 )z — a log 0 < 2 < 1 . 


To find h(z) we let h( oo) 0 (d(oo) ^ 1) and we recall (2.161. Note that <j>%(z) + (j>_>(z) = 0 for 0 < z < 1 

(2.25) 


and (j>ty{z) — <f>^,(z) = 27ri for z > 1. Following [2T], we claim that 


h(z) = — ^alog^ z + (a + 1 )z + ^cnri — a (/>_► (z) — (z(z — 1 )) 


1/2 


is an appropriate choice. First, we have that for z > 1, 

h + {z) = -^cdog z + (a + 1 )z + ^mri - a^4>%(z) ~ V z ( z ~ 1 )> 

h~(z ) = —-a(logz + 27ri) + (a + 1 )z + - a 7 ri — a-((j)t ¥ (z) — 27ri) — \J z(z — 1) = h + (z). 


Therefore h is analytic on C \ [0,1]. So for 0 < 2 < 1 

h + (z) = --alog 2 + (a + 1 )z + -mri - a-</>+(z) - iy/z(l - z), 

h~(z) = —-a(logz + 27ri) + (a + 1 )z + -a 7 ri ~ ot^4>Z,(z) + iy/z(l — z), 
h + (z) + h~(z) = (2a + 2 )z — alogz. 


Finally, lim z _ ) . 00 h(z) 
and then 


ia(21og2+l) + i follows from (2.16). ThenH(z) = e h ^ a3 , D 0 


._ e (alog2+i(a+l))o- 3 


Soo(z) = D~ 1 U(z)D(z). 

We note that we can also write, using (2.13), D(z ) = e h ^ cr3 where 


h(z) = -^alog^ z + (a + l)z - a^4>^(z) - (z(z - 1 )) 


1/2 


(2.26) 


3 The asymptotics of critically-scaled Laguerre polynomials 


From the above calculations, one may guess that S(z) ~ Soo(z) in some sense as N —>• oo. However, this 
is not justified because the convergence of the jump matrix of S to the jump matrix of Soo is not uniform. 
Actually, because of the singularity behavior of S 00 (z) at z = 0, we’ll see that S(z) 76 S OQ (z). We now 
develop local parametrices to solve the Riemann---Hilbert problem for S locally near z = 0 and z = 1. This 
requires the construction of the so-called Airy and Bessel parametrices. 


Remark 3.1. Unless otherwise noted we use the convention 


' M u 

Mi 2 

0(a) = 

' 0(a)M 11 

0(a)M 12 ' 

Ad 21 

M 22 

0(a)M 21 

0(a)Ad 22 


Specifically, Of) should be treated as a scalar, but it can be a different function in each component. 


16 























Figure 3: Dividing the complex plane to define the Airy parametrix with the contour EAi = 71 U 72 U 73 U74. 
The contour 72 makes the angle 7 r /3 with 73 . 


3.1 The classical Airy parametrix. 

First, we divide the complex plane for variable £ into sectors, see Figure[3j Let co = e 27r1 / 3 and Ai(£) denote 
the Airy function. Define 


Ai(£) Ai( w 2 0 
Ai'(0 w 2 Ai'(w 2 <e) 


co 


- 0 - 3/4 

5 


if £ G 1 , 


Pm(0 = 


Ai(0 

Ai'(0 


—co 2 Ai(w£) 
- Ai'(w£) 


CO 


-<t 3 /4 

1 


Ai(0 Ai( w 2 £) 
Ai'(f) w 2 Ai \co 2 £) 


1 0 
-1 1 ’ 


if i g iv, 


if £ G II, 


Ai(0 

Ai'(0 


—co 2 Ai(w£) 
- Ai'(o;0 


w -^3/4 


1 0 
1 1 ’ 


if £ G III. 


(3.1) 


From the asymptotic calculations in Appendix |B.1| it follows 
Hilbert problem: 


that Pa; solves the following Riemann- 
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Riemann Hilbert Problem 3.1. 


PW = Pm(0Jm(0, Z e s Ai , 


Jm(Z) = 


1 l 
o 1 

1 0 
1 1 

0 1 
-1 0 


ifZ e 71 , 

*/£e 72 U 74 , , 
ifZ e 73- 


T A i(£) = 


_ J _ e - 1 / 4 e -|? 3/2 ^£- i / 4 e §€ 3/2 

_i_tl/4„-fj 3/2 <W 4 

2V5F“5 


2V5F 


^ 1 / 4 e 3? 


2t3/2 


w -<W4 ^ + 0(C 3/a )) ■ 


We will come back to this and use it heavily later. We also rewrite the the asymptotics in a more 
convenient form: 


Pm(0 

Em® 



EM®e~ 


2 t3/2 
3^ 


0"3 


Ct2 1 / 4 (1 -)- 3/2)) w 1 / 2 (l + 0(^- 3/2 )), 

_ -W- 1 /4(l + C>(^-3/2^ w 1 /2(1 + 0(f-3/2)) 


3.2 The classical Bessel parametrix. 

The Airy parametrix has the characteristic that four contours exit from the origin in the £ plane. This is 
the case, locally, near z = 1 in our Riemann-Hilbert problem for S. Near 2 = 0 it, the contours look more 
like those in Figure [4] Define 


I a (2£ 1/2 ) - K ct (2£ 1 / 2 ) 

2m£ 1 / 2 I' Q (2£ 1 /2 ) _2£i/2 K ' Q (2£ 1 /2 ) 


TBessel(£) — < 


iHW(2(-£) 1 /2) i H i 2 )(2(-£)V2) 

7r£ 1 /2 H ^ 1 ) , (2(—£) 1 / 2 ) 7 r£ 1 /2 H W' (2 (-£) 1 /2) 


e 4a7ri(T3 


if £ G I U IV, 


if £ G II, 


|Hi 2 )(2(-£) 1 /2) -I H W(2(-£)V2) 

—7r£ 1 / 2 H^ 2 - >, (2(—£) 1 / 2 ) 7r£ 1 / 2 H< K 1 ) , ( 2(-£) 1 / 2 ) 


e"5 a7rilT3 , if £ e HI. 


(3.2) 


Here l a , K a , and are the modified Bessel and Hankel functions [2U]. 
and the jump condition established in [IS] we have the following: 


From the calculations in (B.2) 
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Figure 4: Dividing the complex plane to define the Bessel parametrix with the contour Seessei = 
Here the contour /3i makes an angle of 7r/3 with /? 2 - 


Riemann Hilbert Problem 3.2. When a is given by (1.1 1, the function Pgessei solves: 


p+ 

Bessel 


(f)~ ^Bessel (0 ^Bessel (^)> f € ^Bessel, 

1 0 

^a-Tri ^ 


^Bessel (£) — < 


Bessel 


(0 = 


o 1 

-1 o 

1 o 

e“ a7ri 1 


iff G /3i, 

iff, € /?2, 


iff, G ft, 


O 


C> 


kl “ /2 k |- a/2 

i^i a/2 ir a/2 


ir a/2 

\ f \~ a/2 


ir a/2 

ir a/2 


if I arg£| < 2 tt/3, 


if 2tt/ 3 < | arg£| < tt, 


os £ —> 0, 


PBesse\(M 2 = 


1 I 1 \ 4 / 2 c-l/4.2Md /2 P -cr 1/2 i / 1 j 1 / 2 t-l/4„-2Md /2 r 'Cr 1/2 

2 It tM ) ? C e 2 VttM/ S e e 

i (7 tM) 1/2 ^ 1 / 4 e 2M S 1/2 e - c ^ _1/2 (nM) 1/2 e / 4 e- 2 ^« 1/2 e ^- 1/2 


x (/ + C>(a -i )), Af-»- oo, |^| >1/0 0. 
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As in the case of the Airy parametrix, we rewrite this in a more convenient form: 


P BesS el(M 2 £) 

pBessel(Af 2 £) 


(7rM)-i ff3 r iCT3 ^Be SS ei(M 2 0e 2M?1/2CT3 e- c « V2ff3 , 

'la+0(0) i(i+ 0 (O)i 

i(l + 0(«- 1 )) l + L^a- 1 ) _ ■ 


These asymptotics apply for all £ with | arg£| < n. Furthermore, the asymptotics remain valid up to the 
boundary, arg£ = ±7r. 


3.3 Mapping the Airy parametrix. 


As it stands, the function P A i solves a Riemann-Hilbert problem that resembles the jumps of S near z = 1 
but not exactly. We perform a change of variables and pre-multiply by an analytic matrix function to make 
this match exact. There is an additional constraint. We want the resulting local solution to also match with 
Soo in an appropriate manner. Consider for z > 1 

<j) 4 -(z) = 2 J (s — l)^ 2 (s)jl 1 ^ 2 ds. 

Expanding (s)jl 1 ' /2 = 1 + G(s ) in a power series about s = 1, it follows that 

<M*0 = - 1) 3/2 (1 + (z- l)G(z)), 


where G{z) has a convergent power series in s — 1 with real coefficients. Define the function 

/«-(*) = 2 2/3 (* - 1)(1 + - l)G{z)f/\ 


(3.3) 


which is analytic function in a neighborhood of 2 = 1. Here (-) 2 ^ 3 denotes the principal branch. We then 
have 


/<-(*) = 



2/3 


in the sense that (f^(z)) 3 ^ 2 = § 4><~(z) for 2 close to one. We establish the following facts about /<_: 

• /<_( 1) = 0 and /(_(!) > 0. 

This follows from (3.3). This shows that, in particular, /<_ is one-to-one (conformal) near 2 = 1. 


• f<_(B(l, 5) fl R) C R. for sufficiently small S > 0. 
This follows because G(z ) is real for 2 real. 


• 5) n C 1 * 1 ) C C ± for sufficiently small <5 > 0. 

First, let S' be sufficiently small so that /<_ is one-to-one on B( 1, S'). Let L > 0 be the largest value such 
that (— L, L) C <5 7 ]). Now, let <5 < S' be sufficiently small to ensure that f<-(B( 1, S)) C B{ 0, L). 

Assume that Im/<_(a) > 0 and Im/<_(&) < 0 for a, b £ B(1,S) fl C + . Therefore, on the line that 
connects a to b there must be a value z* such that f(z*) £ (—L, L) and this contradicts that /<_ 
is one-to-one. Then (1, <5) f~l C + ) must be mapped into either the upper- or lower-half planes. 

Considering /<_(1 + ie) = 2ie(l + 0(e)) 3 / 2 we see that Im/<_(1 + ie) > 0 for sufficiently small e and the 
claim follows. 


Consider 


S<_(«) := MAi{z)PAi{M 2 ' 3 f^{z))e^ z)a3 e N<t> ^ z)cT \ 


M Ai (z) := 2y/^D 00 1 J\f(z){'ip<-(z)) a3 


,-i/4 


,!/2 


wn - 1 / 4 u, 1 / 2 


(M 2 / 3 /^(2))* CT3 , 


M := N + -(a + 1). 
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K s 



B(1,S) 


Figure 5: The conformal map £ = f<-(z) applied to B{ 1, 8). We choose T^ so that it coincides with TAi in 
C + and so that it coincides with TAi in C . 



K S (~l 51 Ai 



r Ai 


Figure 6: The pullback of Ks H Eai to create the contour TAi- 


Let 8 > 0 be sufficiently small so that /<_ is one-to-one, analytic and maps C + into C + when restricted to 
5(1, <5). Let Ks = /<_(5(1,5). It follows that I\s is an open neighborhood of the origin. Define a contour 
r A i := fZ- l i K & n s Ai)- See Figures I and § for a graphical representation of this procedure and how this 
affects the precise definition of F| and inside the ball 5(1, 6). 

We prove the following lemma in Appendix |B.3| 

Lemma 3.1. The functions S^-(z) and Mp^{z) have the following properties: 

• S<-(z) is analytic in B(l,5) \ Tai and is continuous up to Tai- S<-(z) has the same jumps as S(z ) in 
a neighborhood of z = 1. 

• For \z — 1| = 5 

S^(z) = D- 1 E M (z)D 00 S 00 {z), (3.4) 

where Ea_\{z) = I + 0(M _1 ) as M = N + |(a + 1) — > oo. 

• Mai(z) is analytic in a neighborhood of z = 1. 

3.4 Mapping the Bessel parametrix. 

We perform the same steps for -Peessei so that it solves the Riemann-Hilbert problem for S(z) near z = 0. 
As before, we would want the resulting local solution to match with S ^ but we will see, this is impossible 
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and this complication will modify the entire discussion that follows. 
Consider for z < 0 


</>->■ (z) = 2 [ (s)^ 1/2 (s - l)^ 2 ds = 2 [ — \/l — sda. 

Jo Jo V~ s 

Now because y/l — s has a convergent Taylor series about s = 0, with real coefficients, we find 

</Mz) = -4 a/R(1 + zT{z)), 

for an analytic function T(z) whose Taylor series has real coefficients. Define 

f^(z):=~4z(l + zT(z)) 2 . (3.5) 

It is clear that this function is analytic in a neighborhood of z = 0. We further note that for z qL [0,oo 
sufficiently small 

C/U(z)) 1/2 = -i<Mz), 

because f->(z) and —<j>^(z) are both positive for z < 0. 

We establish the following facts concerning /_>: 

• /_>.(0) = 0 and /(_(0) < 0. 

This follows directly from the dehnition of /_>. 

• /_>(!?(0, 5) fi K) C K. for sufficiently small <5 > 0. 

This follows from the fact that the Taylor series of \/l — 3 has real coefficients. 

• /_>(!?(0, 6 ) n C 1 * 1 ) C C T for sufficiently small <5 > 0. 

The same argument used to show that f < -(B(l, S) D C^) C C ± can be applied to —/_>.(z) and this 
claim follows. 


Consider 

PBe S s e i(M 2 U(z)) = (TrM)-^(f^(z))-^E Bessel (M 2 f^(z))e- M ^ z ^e- 2c /^ z l 
Note that <(>_>. ( 2 ) only vanishes at 2 = 0 from ( |2.22 1. Consider the function 

S^(Z) = MBe SS elWPBe SS el(M 2 A(^))e^ (z)CT3 e JV ^W CT3 , 

i i l " 1 

M Bessel (z) -.= D^mzW^z))-** \ l (M 2 U(z))^n^, 

w(z) := (2a + 2 )z — a log_^ z + (a + l)7ri. 


(3.6) 

(3.7) 

(3.8) 


In the same way as before, let S > 0 be sufficiently small so that /_> is one-to-one, analytic and maps C + 
into C“ when restricted to B(0,S). Let Lg = f^(B(0,6). It follows that Lg is an open neighborhood of 
the origin. Define a contour Vessel := ffs(Lg fl Vessel)- We reverse the orientation of all contours. See 
Figures [7] and [8] for a graphical representation of this procedure and how this affects the precise definition of 
T r and T|. inside the ball 13(0, S). 

We prove the following lemma in Appendix |B.4| 

Lemma 3.2. The functions S^(z) and M Besse \(z) have the following properties: 

• S-y(z) is analytic in B(0,S) \TB eS sei and is continuous up to Vessel \ {0}. S i -(z) has the same jumps 
and singularity behavior as S(z) at z = 0, that is 

— S^(z) =0(1) as z 0 from outside the lens and 
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Figure 7: The conformal map £ = f^(z) applied to F?(0,(5). We choose T^ so that it coincides with TRessei 
in C+ and Fj, so that it coincides with TBessei in C . 



Lfi 0 SB esse i i fiessel 

Figure 8: The pullback of L$ n Seessei to create the contour FBessei- 
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•SUO) 


i o 

-j- e 2 N<j>^(z)+w{z) ^ 


= 0(1) as z —>■ 0 inside the lens. The (+) sign is taken for z in 


the region enclosed by [0,1] and T-j- and the (—) sign is taken in the region enclosed by [0,1] and 

n- 


(3.9) 


• For \z\ = 6 

S->(z) = ^Bessel (z)D 00 S 00 (z)e 2c ^ z)a3 , 

where ^Bessei(-z) = I + 0(a _1 ) as a —)• oo. 

• MBessei(~) is analytic in a neighborhood of z = 0. 

To account for the fact that (3.9) has a factor of e 2 c/^(z)cr 3 we cons i(jer an additional Riemann-Hilbert 
problem. 

Riemann-Hilbert Problem 3.3. The function A(z) satisfies the following properties: 

1. A(z) is analytic on C \ dB( 0, 6) for 0 < <5 < 1 and continuous up to dB( 0, S), 

2. A+(z) = A-(z)S oc (z)e- 2c /^^ a3 S- 1 (z) for \z\ = 6, and 

3. A(z) = I + 0(z ~ l ) as z —> oo. 

Here dB(0,S) has counter-clockwise orientation. 

First, we look closer at the jump matrix 

S 0 o(«)e“ 2 c/ ^ ( * )os S , - 1 («) = D^M(z)D(z)e- 2c /^ z ^D- 1 (z)M~ 1 (z)D 00 

= D^ 1 Af(z)e- 2c ^ z) ^Af- 1 (z)D 00 , 

because D(z) is diagonal. Now, define A 00 (z) := D 00 A(z)D^f and it follows that A 0 0 has the jump condition 

A+iz) = A- 00 (z)N(z)e- 2c '^^U~ 1 (z). 

Note that A^ is independent of N. This jump matrix is analytic across (0,1): Indeed, for 0 < z < 1 


Af + (z)e- 2c/ ^ {z '> IT3 Af^\z) = Af-(z) 


0 1 

-1 0 


a 2c/0_> (0)0-3 


0 -1 

1 0 


a r\z) 


= N-(z)e~ 2c/(l, ^ z)cT3 Nz 1 (z). 


Following standard theory, the solution is unique if it exists. We now construct an integral representation 
for A ca (z) assuming it exists. Once we have the representation it is a rather simple matter to check that 
the formula gives a bonafide solution. First consider, B(z) := A 00 (z)N(z). Then B(z) is analytic in C\r^, 
T s = dB(Q, 6) U [0,1] and it has the following properties 


B+(z) = B-(z)e- 2c /^^ a 3 , \z\ = 6, 


B+(z ) = B~(z) 
B( 00 ) = I. 


0 < z < 1, 


Let Bi(z) = / 2 (z)] be the first row of B(z). It follows that / 2 (z) = 0(z - 1 ) as z —> 00 and we are led 

to consider g(z) = (z(z — l)) 1 ' 2 f 2 (z)fi(z). Then for some £ £ C to be determined below 

g + (z) = g~(z), \z\ = 6 or 0 < 2 < 1, g( 00 ) = £. 
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Because B{z) should have at most fourth-root singularities at z = 0,1, we require g(z) to be bounded at 
z = 0,1. Thus g(z) is entire, g{z) = £ and hence 


h(z) = 


A(z){z(z-i)y/*- 
We are led to the following Riemann-Hilbert problem for fi(z): 

A + W=/fWe- 2c/ ^ W , \z\=S, 

fi(z)fr(z) = -£(z(z-l))Z 1/2 , 0 < z < 1, (3.10) 

fi (oo) = 1. 

Now, using the principal branch of the logarithm, consider the function r(z) = (log /i(z))(z(z — l)) -1 / 2 : 
r + (z) — r~(z) = —2c/(j)^{z){z(z — l)) -1 / 2 , \z\ = 6 , 
r + (z) - r~(z) = (z(z - 1))+ 1/2 log (~£{z(z - 1))~_ 1/2 ) , 

= (z{z - 1)); 1/2 (log(-tf) - i log(z(l - z))) 0 < z < 1, 
r{z) = 0(z~ 2 ), z — y oo. 

We choose £ to enforce the last condition. The function r(z) is the sum of a function that satisfies the first 
jump condition and a function that satisfies the second condition. We first claim that 


~(z) := - log 2 + 


2z — 1 


(z(z-D)^J’ 

is analytic in C \ [0,1]. This is true as the argument of the logarithm is never negative. Indeed, if 
, z-1/2 


(z(z- l)) 1 / 2 


= -y 2 and so z 2 - z - b = 0, b = 1/4((1 + 7 2 /2) 2 - 1) > 0, 


1 ± Vl + b 

or z = --- € 


But because (z(z — l)) 1 / 2 > 0 for z > 1 and (z(z — l)) 1 / 2 < 0 for z < 0 the claim follows. It also follows 
that r + (z ) + r~(z) = log(z(l — z)) and hence 




2c/</>_,(s) 


{|s|=(5} {S(S~ 1 ))V 2 (S-Z) 


ds+ (log(^ M )_l. w ) (2(2 _ 1)) - 1/2 (3U) 


Then r(z) = r\/z + 0(z 2 ) as z —> oo where 


1 


r i = 


2rri J { \ s \= S } (s(s - l)) 1 / 2 
We must have rq = 0, so we choose £ = £{c) by 


2c/^(») ds+ log(-ig(c)) +log2 


£(c) = iexp [ — c—. [ 

V nlJ { 


2/<M* 


{H = 4} ( S ( S -l))l/2 


ds — 2 log 2 


Because </>_>.(z) = </>~y(z) for z £ [0, oo), the exponent is real, implying that l € iM + . Then 


(z(z- l)) 1/2 r(z;£(c),c) = -4 [ 

7T1 Js 


d» + ilog(-i« C ), + ilog(2 + 
{M=«} (s(s- l)) 1/2 (s-z) 2 2 


2z — 1 

wz-D) 1 ^;’ 

(3.12) 
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Figure 9: The full definition of the contours Tf and F|. The definition inside the shaded regions is given 
by TAi and Tgessei and the remainder is a straight line connecting the two contours. Note that the contours 
are continuous but not necessarily smooth. 


fi(z) = exp ((z(z - l)) 1/2 r(z;e(c),c)) and f 2 (z) = £(c)(z(z - 1)) 1/2 exp(-(z(z - 1 )) 1/2 r(z;£(c),c)). It is 
easy to see that both fi and f 2 have at worst fourth-root singularities at 2 = 0,1. 

To determine the second row B 2 (z) = [h\(z), h 2 (z)\ of B(z) we use the same procedure. It follows that 
hi(z)h 2 (z)(z(z — l)) 1 / 2 = k is constant. We obtain a Riemann-Hilbert problem for h 2 (z) 

ht(z) = h 2 (z)e 2c ^ z \ \z\=5, 
hi (z)h ,2 (z) = k(z(z — 1))_ 1 ^ 2 , 0 < z < 1, 
h 2 (o o) = 1. 


If one sends c —> — c and k —> — £ then we have (|3.10|). So, h 2 (z) = exp((z (2 — l)) l / 2 r{z] — £{— c), —c)) and 
h^z) = -£{-c){z(z - 1))" 1/2 exp(-(z(z - 1 )) 1/2 r(z; -£(-c), -c)). Finally, we have 


Aoo(z) - 


e (z(z-i))^ 2r ( z; e(c),c) 

-£(-c)e-^ z - 1 '>'> 1/2r( ~ z '-^- c ^~ c \z(z - 1)) _1/2 


£{c)e-W z - l » r( 2 ;t( c),c)( z ( z _ l))-l/2 

p ( 2 ( 2 -l)) 1 '' 2 r(z;-n-c) ,-c) 


M~\z), 


and using (3.12) we see that det^4oo(z) = 1. 


3.5 The final deformation 

Before we discuss the final deformation, we must specify T-f and F|. Based on the discussion in the previous 
section we have defined both contours TBessei and Tai inside the balls i3(0, <5) and B(1,S), respectively. Let 
r r and F| to be the contours obtained by connecting TBessei and Tai with straight lines as in Figure[9] For 
convenience we also define = T-j- \ (B( 0, S) U B( 1, 6)) and similarly T^ = Tj. \ (13(0, S) U B(l, S)). 

Define 

( D 00 S(z)S^ 1 (z)A~ 1 (z)D if Z G B( 0, S) \ r B essel, 

£(z) = \ D 00 S(z)S- 1 (z)A-\z)D^\ if Z€B(1 , 6) \Tm, _ 

{ D 00 S(z)S~ 1 (z)A- 1 (z)D-\ if * € C \ (13(0,5) U 13(1, 6) U (0, oo) U T t U T x ). 

We look at the jump of £ on <913(0, 5): 

£+(z) = D 00 S(z)S~ 1 (z)A-\z)D- 1 

= D 00 S(z)S- 1 (z)A+ l (z)D- 1 (D 00 S(z)S- 1 (z)AZ 1 (z)D- 1 )- 1 D 00 S(z)S- 1 (z)AZ 1 (z)D 
= £-(z)(D 00 S(z)S- 1 (z)Az\z)D^)- 1 D 00 S(z)S- 1 (z)Al\z)D^ 

= £-(z)D 00 A_(z)S 0C (z)SZ>\z)A+ 1 (z)DZ 0 1 . 
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From the jumps of A(z) we find 


£ + (z) = £-(z)D 00 A+(z)S 00 (z)e 2c / <l, ^ a3 S- 1 (z)A^(z)D^ 1 . 


Using (3.91 


£ + (z) = £-(z)D 00 A + (z)D^E^ ssel (z)D oc A^(z)D~ 1 = £~ (z)A+(z)E^ ssel (z)(A+ )~\z). 


(3.13) 


For z € 313( 1, 6) 

£ + (z ) = D oc S(z)S- 1 (z)A-\z)D- 1 

= D 00 S(z)S- 1 (z)A-\z)D- 1 (D 00 S(z)S-\z)A- 1 (z)D- l )- 1 D 00 S(z)S-\z)AzHz)D 
= f-(2)( J D 00 5(^) ( 5- 1 (z)^- 1 (z)D- 1 )- 1 ^ 00 5( 2 )5- 1 (^)A- 1 (z)D- 1 
= f-(2)D 00 A( 2 )5 00 (z)5- 1 (2)A- 1 (z)i?- 1 . 


Using (3.41 


£ + (z) = £-(z)D 00 A(z)D- 1 E^(z)D 00 A~ 1 (z)D- 1 = £~ (z)A+(z)E^(z)(At 0 )~ 1 (z). 


(3.14) 


From here on, we let 0 < <5 < 1/2 be sufficiently small so that that /«_ and /_> are one-to-one on 13(1, 8) 
and 13(0, 8), respectively. We establish the following: 

Riemann Hilbert Problem 3.4. The function £(z) satisfies the following properties: 

1. £{z) is analytic on C \ T,r for T,r = F) U T) U 313(0, <5) U 313(0, 5) U (1 + <5, oo) and continuous up to 
the contour T,r. 

2. The jump conditions for £(z) are given by 


£+(z) 

Js{z) 


£ (z)Js{z), ZG'Er, 

I I + «/ z eE fi \3B(0,i), 

\ 1 + C>(a- 1 ), if zG 35(0, <5), 


where the error terms tend to zero in L 2 n L°° as N —» oo. 


3. £{z) = I + 0(z ) as z —► oo. 


Using (3.14) and (3.13) and Lemmas 3.2 and 3.1 it is clear that the jump matrix for £ has the correct 
behavior on 313(0, 6) and 313(1, S). On Ti, r( and (0, oo) the jump matrix for £{z) is given by 

A 00 (z)Af(z)D(z)Js(z)D~ 1 (z)Af~ 1 (z)A^(z). 


Because A 00 (z),A 0 f(z),J\f(z) and Af x {z) are all uniformly bounded on the contours being considered we 
only need to consider D(z) Js(z)D^ 1 (z). For z > 1 + <5 


D(z)J s {z)D 1 (z) 


^ g—2 N(f>~^ (z)-\-a7r—a(j>~^ (z) —2(z{z— l)) 1 / 2 

0 1 


We have Re0+ (z) > C > 0 for 2 > 1 + <5. Let R > 0 be so that </>+( 2 ) > n for 2 > R. Then on [1 + 5, R] 
we have 


< e 


-2NC 


0 in L 2 nL°°, 


exp (—2Ncj>f^(z) + an — a<j)f f {z) — 2 ( 2(2 — 1)) 1,/2 ^ 
for any e > 0. On [R, 00 ) we have 

exp 2N(j)fz{z) + an — acj)fj.(z) — 2 ( 2(2 — l)) 1 ^ < exp (—21V Re <f>t^.(z)) —> 0 in L 2 0 L°°. 
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On U T[ we have 


D(z)Js(z)D~ 1 (z) 


1 0 ' 

e 2iV0+(z) — a7r+a0+ (z)+2(z(z—l))]^ 2 ^ 


From Lemma 2.2 there exists a (possibly new) constant C > 0 so that Re<jf>j),(z) < —C on Ti U T(. Let 


m = max r ^ ur ^( 2 | 2(2 — 1)| 1//2 ) then it follows that 


exp (2N(/)ty{z) — an + a<f>%(z) + 2 (z(z — l))^ 2 ^ <e ( 2N + a ) c —> Q in L 2 n L 


Finally, we have to establish that £(z) is analytic in B(0,8), B(l,6) and in a neighborhood of each point in 
(<5,1 — <5). We compute the jumps for 2 £ Vessel • First 

£ + (z) = D 00 S + (z)(S+)- 1 (z)A~ 1 (z)D^. 

Then from Lemma [T2] 5 ,+ (z)(5 l j > )~ 1 (g) = S~(z)(S^)~ 1 (z) and this is an analytic function in B{ 0, 5). Thus 
£ + (z) = £~{z). Similar calculations follow for 13(1,5) and (5 ,1 — 5). 

From classical theory (3] (see also [2B]) it follows that 

£(z) = I +Cy. r u n {z), u n = 0(a _1 ) -t 0 in L 2 (Y R ). 


It then directly follows that £{z) —> /, £'{z) —> 0 uniformly for 2 bounded away from £#. We have thus 
proved the following theorem concerning the asymptotics of Laguerre polynomials. 

Theorem 3.1. There exists <5* > 0 such that for any 8 £ (0,5*), the upper limit, Y + (x), x £ (0, oo) takes 
the following form: 

(a) In the region 0 < x < 8, 


Y + (x) = e~^^D^ 1 £{x)A 00 {x)D 00 SX{x) 

(b) In the region 1 — 8 < x < 1, 

Y+(x) = e~^ (Ts D^ 1 £(x)A 00 (x)D 00 SUx) 


1 

g2N4>~f(x)+w(x) 


1 

g2N(p'f(x)-\-w(x) 


0 

1 


0 

1 


e (Ng + (x)+±e N )a 3 


e (Ng + (x)+ i£jv)cr 3 


(c) In the region 1 < x < 1 + 8, 

Y+(x) = e-^ lNa3 D^£{x)A 00 {x)D 00 St{x)e {N9+{x)+¥N)a3 . 

(d) In the region 1 + 6 < x < oo, 

Y+(x) = e-i^D^£(x)A 00 {x)D 00 S+{x)e^ Ns +^ + i eN ^ 3 . 
Further, as N — ► oo, £{x) — > I, £'(x) —)• 0 uniformly in each of the regions listed above. 


4 The extreme eigenvalues and the condition number 


In this section, we prove Theorem 0 Theorem |1.2| and Theorem |1.3[ The proof of Theorem 0 and 
Theorem |1.2| relies on the asymptotic results of Section [3j and some basic operator theory. Since the kernel 
ICn is related to Y through equation ( |2.7[ ), we apply Theorem 3.1 to show that after suitable rescaling, the 
kernels ICn converge to the Airy kernel JCai at both the soft and hard edge. Standard results on operator 
theory are then applied to establish convergence of the associated Fredholm determinants. The proof of 
Theorem |1.3| follows from these results, though an additional lemma is necessary to account for the fact that 
the largest and smallest eigenvalues are not independent. 


28 












4.1 Uniform convergence of /Cy at the hard edge 

We first consider the smallest eigenvalue. Let <5* be as in Theorem |3.1[ fix S £ (0, <5*) and consider x in the 


interval (0, <J]. We use equations (3.6)—(3.7) relating S to the Bessel kernel and Theorem |3.1[ a) to obtain 
1+(U = e~^ Na3 D^ D 1 £(x)A 00 (x)D 00 M B e SS e\{x)P B e SS e\(M 2 


1 

e -(a+l)7ri 


e cV</>+ (x) +Ng + (x)+ i l N )<r 3 e § w + (x)ct 3 


It follows from the definition of MBessei(z) that MBessei(^) : = D<x>M B essei{z)M 2 CT 3 j s an analytic function 
in 13(0, (5) that is independent of N (also of M and a). It is convenient to introduce the following functions: 


V{z) = w3(z)[ 0 1 ] 


1 0 

Q:7ri ^ 


[PBessel(M 2 A(0))_]- 1 M-5-3 ) 

(4.1) 


W(z) = M^P Bessel (M 2 f^(z))_ 


1 

3 —(a+l)7ri 


3 ( N <t>X (•*)+ Ng+ (z)+ \t N )cr 3 e I Jl>+ (z)cr 3 


W3(z). (4.2) 


A straightforward calculation then yields 


w3(x)w3(y)[ 0 1 ] Y^(x)Y + (y) [ 1 0] = V(x)M^ sael (x)A^(x)£- 1 (x)£(y)A 00 (y)M Bessel (y)W(y). 


Thus, we may define the analytic function B(y) := £{y)M Besse i(y) and rewrite equation ( 2.7) in the form 

K- , n 1 V{x)B-\x)B{y)W{y) 

A N{x,y = - — -. 

27ri x — y 

The study of the asymptotics of K N has now been reduced to the asymptotics of V and W. We further 




simplify V and W using the definition and properties of the Bessel parametrix. Comparing (3.2) and ( |4.1| ), 
we see that we must consider the Bessel parametrix PBessei(C) with argument £ = M 2 /_ > (x) with x £ (0,1). 
The relevant regime here for the Bessel parametrix is (III). Indeed, for x £ (0,1), </> + (cc) is purely imaginary 
with a positive imaginary part so that f->(x) < 0. Furthermore, because /_> maps the upper-plane to the 
lower-half plane (at least locally) (negative imaginary part) and (—C) 1 ^ 2 = — ^M<j>t ¥ {x) 

(positive real part). Recall also that Pscssei easily computed since pBessei has determinant one. From 
( |3.2| ) (III), we obtain 

-%M<l>+(x)HW\-iM<j>+(x)) \ R^i-iMcfr+ix)) 

-f Mcp+ (x)H^'(-i^i(U) 5 H^(-iM0+(*)) 

= e ( N + 2 ') 7ri M 1//2 [ —n(j)i y (x)J ol '(—iM(t>t > (x)) J a (—i M(f>t+{x)) ]. 

In the second line we have used the Bessel function identity J a (a;) = |H^(a:) + \ H^ 2 \x). A similar 
calculation for W yields 

i H«(-i M<t>+(x)) -|H«(-i M^+(x)) 

§Af0+ (x)U^\-iM^+ (a:)) -^M<j>+(x)B^\-iM<l>+(x)) 

e ( N +h h'M 1 / 2 


V{x) = e {N+ s )wi [11] 


(4.3) 


W( x) = M? c 


1- 

1 

-1 

-1 


a (Af+i)7ri 


Ja(-i M<j>%.{x)) 
’n(j)t f {x)3 <x '{-\M(j)X (x)) 


(4.4) 


The above formulas hold in the region x £ (0,1). It is convenient to extend the range of definition of V and 
W by adopting the convention V(x) = W(x) = 0 for x < 0. 

We now apply classical asymptotic formulas for the Bessel function J Q to obtain quantitative decay 
estimates on V and W. These will imply convergence results for /Cy. More precisely, define the rescaled 
variables and the rescaled kernel 


c 2 2 2/3 


or 


, 8/3 


1 + X 


or 


2^/3N 


£n(x, y)dy := K. N (x, y)dy. 


We then have the following convergence result. 
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Proposition 4.1. As N — > oo the rescaled kernels converge pointwise, 

c- , \ Ai(— y) Ai'(— x) — Af(— y) Ai(— x) . . 

JC N (x,y)^^^ -i- L - v y> K - (x,y)£ K 2 , 

x-y 

and the convergence is uniform for ( x,y) in any compact subset of (—oo,L] 2 for any L £ K. If x = y 
then the limit is determined by continuity. Further, there exists a positive, piecewsie- continuous function 
G : (— oo,L) 2 —>• (0, oo), such that 


\^N(x,y)\ < G(x,y), 


G(x, y)dx d y < oo, 


G(x, x)dx < oo. 


(4.5) 


' —oo J —OO 


The technical lemmas that underly this result are stated below, but proved in Appendix [C] Let W = 
\W\,W 2 } t and V = [Vi, V 2 ]. 


Lemma 4.1. Define 


/ \ 1 e 1 l^+il 372 j/ — 00 < x < — 1, 

9(I) = 1 1 , othcrw.se. ~ 


Assume x £ (— 00 , L], L £ K. Then there exists a constant Cl > 0 and A > 0 such that if a > A then 

\Wi{x)\ = \V 2 {x)\ < C L M 1 ^ 2 a~ 1 ^ 3 g(x), 

\W 2 (x)\ = \V t (x)\ < CLM-^a^gix), 

| W[{x)\ < C l M 1 / 2 a 7/3 , 

\W^{x)\ < C l M - 1 / 2 a 9/3 , 

with the convention that W(x) = V(x) = 0 for x < 0. 

Lemma 4.2. For x in a compact subset of (— 00 , L], L £ R. 

2/3 / 2 \ !/ 3 


V{x) = e (Ar+ 3 )7ri M 1/2 
W{x) = e (Ar+ 3 )7ri M 1/2 

W\x) = e (JV+ s )™M X ' 2 


.a (2 

l M 


(M\~x) + 0(a-^ 3 )) (-J (Ai(-z) + ©(a" 1 / 3 )) 


(4.6) 


/o \ 1/3 

(^J (Ai(—x) + 0{a~ 1 ^ 3 )) 


a (2 


2/3 

M V ol ) 

v 3 / o \ 2 / 3 


(Ai .'(—x) + 0(a 1 / 3 )) 


+0(a~ 1/3 )) 

c z M 


(4.7) 


(4.8) 


where the error terms are uniform in x. 

4.2 Uniform convergence of at the soft edge 

We write the expansion for 1 < z < 1 + <5 

Y+(z) = e-^ NrT3 Df 0 1 £(z)A 00 (z)D 00 M Ai (z)P Ai (M 2 / 3 f^(z))+e { ^ z)+N ^ {z)+N9 + {z ' )+ ^ N)a:> . 


From (2.12) and (2.13) we have 


^w(z) + N(j>^(z) + Ng+{z) + ^i N = ^ w(z ) + N<j>+.(z) + Ng+(z) + - Ni tt 

11 _i 

= ^ vz - 2 a 1 °§ z = 2 {z). 
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We simplify these expressions using 

Y + (z)w 1 3 a3 (z)=e-i e »‘ 7 *D^ 1 £(z)A 00 (z)D 00 M Ai (z)M-i^M^P Ai (M 2 ' 3 f^(z)) + 

and noting that B(z) := £{z)A 00 {z)D ao Mxi{z)M~^ rT3 has no ./V dependence, is analytic on (1 —5,1 + 5) and 
has a constant determinant. In view of l2.7t) we define 


V(z) := [ 0 1 ] P^(M 2 / 3 f^(z))+M-^, 

W(z) := M5 ff3 PAi(M 2/3 /^(z)) [1 0 ] T . 

It follows from [201 (9.2.8)] that det PAi(M 2 / 3 /«_(z)) = (27ri) _1 so that 
V(z) = 27ri [ 0 1 ] 


w 2 Ai(w 2 M 2//3 /^(z)) -Ai (co 2 M 2 / 3 f^{z)) 
— Ai'(M 2 / 3 /<_(+)) Ai(M 2 / 3 /^(z)) 

= 27ri[ -Ai '{M 2 / 3 f^(z)) Ai (M 2 / 3 f^(z)) ]M-3 CT3 , 
Ai(M 2 / 3 / < _(^)) 




W(^) : ^ Ai'(M 2 / 3 f<_(z)) 

By analytic continuation, the same formula holds for 1 — 5 < 2 < 1. And then 

F(a:)B _1 (a;)B(2/)W(j/) 


K-n(x, y) = 


x-y 


for appropriate values of x and y. 

Next, for z > 1 + 5 we examine 

Y+{z)wd a \z) = e~^ Na3 D^£(z)A 00 {z)N{z)D{z)e {N9 ^ z)+ ^ N ^wl a3 {z). 
Define the ^independent function B(z) = £(z)A ao (z)Af{z) and the two functions 

W{z) := D(z)e {N9 + iz)+ ^ N)a3 wl a3 {z) [ 1 0 ] T , 

V(z) := [ 0 1 ]w5 a3 (z)e-V ,s +M + t tN)aa D- 1 (z). 

And then 

k - , ^ V{x)B-\x)B{y)W{y) 

lCN{x,y) = -, 

x-y 

for appropriate values of x and y. Consider the scaling operator 

x 

x = 1 + 


(4.9) 


2 2 / 3 M 2 / 3 ’ 

So that 

M 2 / 3 f^(x) =x + e>(Af 2/3 ), 

uniformly for x in a compact set. Define K,N(x,y) through the equality K.N(x,y)dy = /Cjv(i,y)dy. 

Proposition 4.2. As N —> oo the rescaled kernels converge pointwise, 

- . . Ai(x)Ai'(y) — Ai'(x)Ai(y) , 

/CjvO r,l/)-^—^-—- V ; , (x,y)eR 2 , 

x-y 

and the convergence is unifoi'm for (x,y) in a compact subset of [L, oo) 2 for any Lei. If x = y then 
the limit is determined by continuity. Further, there exists a positive, piecewise-continuous function G : 
(L, oo) 2 -+ (0, oo), such that 


fOO pOO 


\^n(x,v)\ < G(x,y), 


G(x, y)dx dy < oo, 


G(x, x)dx < oo. 


(4.10) 
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As in the previous section, we have some technical lemmas that are used in the proof of this proposition. 
All details can be found in Appendix [ d| Let W = \W\,W 2 \ T and V = [Vi,V 2 ]- 

Lemma 4.3. Define 


g(x) 


(i + M) 1/4 { i” 


if — oo < x < 0, 
otherwise. 


Assume x £ [L, <5M 2 / 3 2 2 / 3 ], L £ R. Then there exists constants Cl > 0 and A > 0 such that if a > A then 

\W 1 (x)\ = ^\V 2 (x)\<C L M 1 / 6 g(x), 

\W 2 {x)\ = ^|Fi(z)| < ClM-^x), 

1^(4)| < C l M 5 / 6 , 

\W 2 (x )| < C l M 1/2 . 

Additionally, for x £ [8M 2 / 3 2 2 / 3 ,oo) 

||lF(x)|| = ±\\nm < e~ 2l/3NM ~ 2/3 * < g(x), 


for sufficiently large N. 


The following follows directly from (4.91. 

Lemma 4.4. For x in a compact subset of [L, oo), L £ R. 

V(x) = 2iri [ -M-V 6 (Ai'(a:) + 0{M~ 2 / 3 )) M 1 / 6 { Ai(ar) + 0{M~ 2 / 3 )) 

[ MV 6 (Ai(x) + C?(M- 2 / 3 )) 

W ^ [ M-V 6 {Ai'(x) +0{M~ 2 / 3 )) 

W '(r) = ? 2 / 3 \ ^ 5/6 ( Ai '(^) + 0(M~ 2 / 3 )) 

[ > [ M 1 l 2 {xAi{x) + 0{M~ 2 / 3 )) 


where the error terms are uniform in x. 


4.3 Proofs of the main theorems 

Our main tools for our proofs are from [22] : 

Theorem 4.1 (Theorem 3.4). The map A H > det(/ + A) is a continuous function on J\ and 
| det(/ + A) - det (I -B)\<\\A- expdlAh + \\BWj. + 1). 

Here J\ is the set of trace class operators with norm 

||A||i = tr VA*A. 

Theorem 4.2 (Theorem 2.20). Suppose A n —► A, |A„| —► \A\ and |A* | —> |A*| all weakly and that ||A„||i —> 
||A||i, then ||A — A„||i —> 0. 

From [3] it is known that 

det(J - /Cjv|z. 2 ((t, s ))), 

gives the probability that are is no eigenvalues in the interval (t,s). Thus taking into account the initial 
scaling by v 

IP (An in ft = det (I K*N | Z. 2 ((0,/:) ) ) ? 

^(^max/^ — t) = det (I Ajv|.L 2 ((f,oo)))' 
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The same statements hold for K.n{x : u) and Ay(:r,y). 

Proof of Theorem \l. 1\ We begin with the observation that 

det(/ - AyI x, 2 (( 0 ,*))) = det(/ - Ay | L 2 ( (_( a/2) 2/3 it ))) = det(/ - Ay l^.^))) 
when we use the convention that Ay(x, y) = 0 if x < 0 or y < 0. Define 

K.Ai{x,y) = —K-m{—x, —y). 


Let/, g G C°°(R) with compact support. Then from Proposition 4.1 \fCN(x,y)f(x)g*(y)\ < G{x,y)\\f{x)\\g(y)\ £ 
L 1 ((—oo, t) 2 ). Thus, by the dominated convergence theorem 


iC N (x,y)f(x)g*(y)dxdy 


fCAi(x, y)f(x)g*(y)dxdy, 


' — oo J —OO 


' —oo J —oo 


as N — > oo. Therefore K,n\l 2 {{- oo,t)) —► ^Ai|i 2 ((-oo.t)) weakly. Additionally, convergence in trace norm 
follows: 


/ Ay ( 2 , x)dx —>■ / Ka\{x, x)dx. 

-oo J —OO 

From Theorems 14.21 and 14.II we have that 

w dn^ det (/-£7v| L 2((°y))) = det(/-/C Ai |L 2 ((o,t))) = F 2 (-t). 
Stated another way, for t G R 


lim P ( A min > 

N—>oo 


to¬ 


ot* 


c 2 2 2 / 3 \ 

) 


= F 2 (-t). 


Also, v = 4 A + o(a *) = or /c so that this result can be simplified. We use the following lemma with 



o 

C V 

Cy = 

, Cy = C, 


or 


c 2 v 2 2 / 3 - 

dy = 

a 8/3 > d * 


Lemma 4.5. Suppose P((Ay — Cy)/dy < t) —► F(t) where F is continuous at t. If (cy — Cy)/dy —t 0 and 
dy/dy -> 1 i/ien P((Ay — cy)/dy < t) —► F(t). 


Proof. First, consider 


Ay — Cy 
dy 


<t = 


Ay ~ Cy dy Cy — Cy dy < 


dy dy 


dy dy 


For e > 0, let A"* > 0 be sufficiently large so that |(cy — cy)/d„| < e and |1 — dy/dy| < e. Then 


lim sup P 

N—>oo \ 


/ X N — Cy 
V dy 
/ Ay — Cy 


< t ) < lim sup P ( Xn , CN 
y-v oo v dy 


< 


lim inf P , 

oo y dy 


< t ) > lim inf P 
N—>oo 


Ay — Cy 
dy 


< 


1 + sign(t)e 
t 

1 — sign(t)e 


<l + e) 

1 - e 

e(l + 0 

1 - e 


= F 

= F 


t _l_ e (! + e ) 

1 + sign(f)e 1 — e 

t e(l + e) 

1 — sign(f)e 1 — e 


Because e is arbitrary, the lemma follows. 


□ 
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To see that our choice of c/v,cjv,<ijv and d n fits the hypotheses of this lemma, note that 

a = V4cN + 0(1) => a 2 =4 cN + 0{VN), 
cv = 4cN + 0(\/N), va~ 2/3 = 0(N 5/6 ), 


and then 


The theorem follows. 


cn - cn 
dN 


cv — a 2 
— C c 2„2 2 / 3 
a 2/3 


dN ot 2 
dN cv 


OiN- 1 / 3 ), 


□ 


Proof of Theorem 1.2: We begin with the observation that 

det(J - /C iV |z J 2 ((i ]00 ))) = det(/ - /C/v|i, 2 ((t,oo)))- 
Following the arguments in the proof of Theorem |1.1| we have sufficient conditions for 


lirn P (A max /^ < i) = hrn det(/ - lC N \L^((t ,<»))) = F 2 (t). 

N-t-OO N-¥ oo 


Written another way, using that v = 4M, 

^max V 


lim P 

N—too 


2 2 / 3 M 2 / 3 


< t = lim P 


^max V 


N— >oo \ r/ 1 / 3 2 2 / 3 


< t = F 2 (t). 


This proves Theorem 1 1.2 


□ 


Before we prove Theorem 1.3 we prove a critical lemma from first principles. 


Lemma 4.6. Assume two sequences of random variables ( X n ) n >a, (Y n ) n >o and two sequences of real num¬ 
bers (a n ) n >o, (b n ) n >o satisfy the following properties: 


• Y n > 0 a.s., a ni b n > 0, 

• Y n = a n + b n Y n so that Y n —> t; in distribution, and 


• a n /b n —► oo and — 1| —► 0 in probability. 

Then 





in distribution. 


Proof. We then claim that for each 1st (where t is a point of continuity for F(t)) 

~\F ~ —i u f) 1 ~ 

Y n CLn H - o n t J 

where F(t) = P(£ <t). To see this, compute 


(4.11) 


Xn < _1 

Y n a■ 


n i 


b n t 


X„ 


< 


a n + b n Y n a n + b n t 


-b l -(X n -l)+X n t<Y n 

Or). 
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For 1 > e > 0, consider 

"a. 


(X n -1 )+X n t<Y n = 


+ 1 


^(X n - 1) + X n t < Y n ,^\x n - 1| > e 

On b n 

^(X n - 1) + X n t < Y n , °^\X n - 1| < e 

On. On. 


It is clear, by the fourth assumption that the first term here vanishes as n —> oo so we concentrate on the 
latter. It is certainly true that for t > 0 and sufficiently large n 

limsupP ( j^-(X n — 1) + X n t < Y n , \ X n — 1| < e J < limsupP (—e + (1 — e)t < Y n 

= l-F(-e + (l-e)t). 

Then we use the fact that P (A (~l B) = P(A) + P(£?) — P(A U B ) to find 

A n = { ^(X n -1) + x n t < y„ j, 

B n = \^\X n -l\ < 

^(X n -l) + X n t<Y n ,^\X n -l\ <e) =P (A n HB n ) 

On 0 n J 

= P - 1) + Xnt < Y n ^j + P (B n ) - ¥{A n U B n ). 

It is also clear by the third assumption that linx^,*-, P {B n ) = lim,,,^;*, P (A n U B n ) = 1. We use the estimate 


Therefore 


lim inf P 

n—>o o 


t~(X n — 1) + X n t < Y n ) > P(e + (1 + e)t < Y n ). 

On 


^(X n ~l)+X n t<Y n ,^\X n ^l\<e) > 1 — F(e + (1 + e)t). 
bn O n 


We have shown that 

/ x 1 

1 — F(e + (1 + e)t) < lim inf P I —(- < 

n—> oo ' ' 


hn Un T b n t 


< limsupP [ —— < 


I n a n -\- b n t 


< 1 — F(—e + (1 — e)t). 


Letting e j 0 demonstrates the claim. For t < 0, this argument can be adapted by replacing (1 ± e) with 

0 ? 0 - 

We now modify things and consider for t being a point of continuity of F(t) 


Xn 

Y n 


< — I 1 - — t 


We note that (for fixed t) 


Define 


- 1 ~~t = 


1 


O'n H - 


b 2 

at 


x n =X n [l- -f t 2 

n z 


-1 


-1 


1-ft 2 I — 1 + E n (t), 


h 2 1 

E n (t) = %t 2 


«r i - %i 2 ' 
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and consider 


^{X n - 1) = - 1) + ^E n (t)X n . 

On On 0 n 

It then follows that \X n — 1| ! ™ b 0 and from this we may apply (4.11) with X n replaced by X n to state 
that for any s where s is a point of continuity of F(s), 


Xu < 1 

Y n Cln 4~ b n S 


l-F(s). 


If we set s = t we find 


X r 


I n @n 


—- < — 1 - — f =p U^< 


x„ 


1 


1// o, n d - b n t 


1 -F(t). 


This proves the lemma. 


Et 


Proof of Theorem\1.3\ The proof of this theorem relies critically on on Lemma |4~6| We note that 


'Enin — C C ( 


2/3 

Y n 

where Yn — > —£goe where P(£goe < t) = F 2 (t). We let 6/v = c(2/a) 2 ^ 3 /v, aw = c/v, Yn = A m i n /z-' and 
Xn = Amax/v. Then for e > 0 


P^|X n -l|>e) =p(c0 2/3 |A m >-l 


- II > e 


A m ax - ^ . eis 2/3 /a\- 2 / 3 


> 


j,1/3 2 2/3 - c2 2/3 


(I) 


< 


cv 


,2/3 


Amax ^ 

1/1/322/3 - c2^/3 


(?) 


-2/3 


For any L > 0 there exists iV* > 0 such that for N > N* 


v m /^Amax - ^ ^ e^ 2/3 fa\~ 2 / 3 \ ( A max — V \ 

*“T Ps ^UJ J sb r£? p («? s£ j- 


because v/a —> oo. A similar estimate follows for 


Amax V 


< -- 




,2/3 


„l/3 2 2/3 - c2 2/3 


(?) 


-2/3 


Then because A is arbitrary we have 


AT—>-oo 


lim P -H|X n -l| >e =0, 


or I s -1 A„ — 1| P -^> b 0. Then Lemma 


4.6 


AT—»-oo 

or because a N 2 b N = c l {2/a) 2/3 v 1 


implies 

bN a N 


Xn _L \ 

lim P | Xk - V*- >-t) = 1 - F 2 (-t), 


' x. 


lim P ^^ < t = F 2 {t) = lim 




N—±oo 


bNO-N 


N — yoo l C l v{2/a) 2 / 3 


< t 


This proves the theorem. 


□ 
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A Motivating numerical calculations 


In this appendix, we discuss the simulations of the halting time T e ,E,N,n that motivated the study of the 
critically-scaled Laguerre Unitary Ensemble. Given M samples of the ensemble E = (F, F) which each 
consist of an N x N matrix A and an iV-dimensional vector b, we can compute M samples drawn from the 
distribution of T £t E,N,n- Recall the definition of the fluctuations 


T e ,E,N,n 


T e ,E,N,n ~ ^[T e ,E,N,n] 
\/Var [T e ^E,N,n\ 


T e ,E,N,n ~ 

&e,E,N,n 


( T e ,E,N,n) 


(A.l) 


where (T £i e,n,ti) and &e,E,N,n represent the sample average and sample standard deviation, respectively, 
taken over the M ~^> 1 samples. We plot the histogram of r et E,N,n in Figure 10 for three choices of F. This 
computation indicates universality for T e E,N,n- 


Notation. When F is a Bernoulli random variable, taking values ±1 with equal probability, we call the 
resulting ensemble the positive definite Bernoulli ensemble (PBE). We also refer to the pair E = (F,F) as 
PBE (or LUE of F ~ X c ). Here F is understood to be uniform on [—1,1], 



Halting Time Fluctuations 


Figure 10: A demonstration of universality for T Et E,N,n when n ~ N+ \ \/ 4c N\ with c = 1. The computation 
is taken over 32, 000 samples drawn from the ensemble E. This figure contains three histograms. One for 
F being a Bernoulli random variable taking the values ±1 with equal probability (PBE), one for F being a 
standard real normal random variable and the last for F being a standard complex normal random variable 
(LUE). 


A.l Ill-conditioned random matrices 


In this section we consider the distribution of T e ^E,N,N^ Le. n = N in the case of LUE and PBE. The 


limiting distribution for the condition number is given in (1.2) for LUE. We plot a simulated histogram for 
the condition number when N = 100 and when N = 196 in Figure [ll(a)| again for LUE. In Figure [Tl(b)| we 
plot the corresponding simulated halting time distribution once again for LUE. The computed moments are 


shown in Figure 11(c) for LUE and PBE and indicate that the fluctuations are not universal (compare with 
Figure 13(c) below). 


A.2 Well-conditioned random matrices 

Here we consider the distribution of T c .e.n. 2 N, Le. n = 2 N in the case of LUE and PBE. It is known that 
the distribution of the condition number has a limit with finite mean. This is demonstrated in Figure [l2(a)| 
for N = 100 and N = 196 for LUE. A simulated histogram for T c> e,n, 2 N is shown in Figure 12(b) for LUE. 


From this plot, it is apparent that the discrete nature of the distribution will persist as N —> oo in agreement 
with the discussion in the introduction. 
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n = N n = N 



LUE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

PBE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

100 

282.442 

1290.84 

-1.75573 

5.43125 

100 

279.86 

2124.82 

-0.544804 

6.54769 

200 

643.592 

2809.2 

-2.47523 

10.6812 

200 

641.935 

7069.64 

-0.0881511 

14.4649 

300 

1018.9 

4220.1 

-2.80277 

14.2713 

300 

1021.09 

14353.2 

1.38426 

36.2729 

400 

1404.26 

4718.35 

-2.12134 

23.1877 

400 

1409.61 

28189.1 

11.4373 

706.288 

500 

1786.33 

6390.55 

-2.80864 

21.4515 

500 

1802.01 

41838. 

5.62383 

136.943 

600 

2173.88 

7846.43 

-2.61028 

20.5573 

600 

2197.1 

55448.6 

4.59607 

95.9881 

700 

2565.14 

9163.6 

-1.93369 

33.4126 

700 

2590.53 

70960.4 

4.02138 

52.7732 


(c) 


Figure 11: Ill-conditioned matrices. Numerical calculations when N = n. All calculations are taken 
over 32,000 samples, (a) Histograms for the simulated condition number for LUE when N = 100,196. (b) 
Histograms for T e} E,N,N for e = 10 -14 for LUE. (c) A table of the computed mean, variance, skewness and 
kurtosis of T Ei e,n,n for both LUE. This table indicates that the kurtosis does not converge as N increases. 
Thus, there is no limiting distribution function for the fluctuations. 
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n = 2N n=2N 




(a) 


(b) 


LUE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

PBE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

100 

73.2159 

1.15328 

0.0759259 

3.0488 

100 

72.2804 

2.17711 

0.0195341 

3.05985 

200 

85.5885 

0.982886 

-0.0105681 

3.03685 

200 

84.9429 

1.88208 

0.00520835 

3.01454 

300 

90.5032 

0.790077 

0.00742551 

3.03203 

300 

90.0251 

1.51354 

0.0114286 

3.031 

400 

93.2354 

0.727263 

0.0189789 

2.96987 

400 

92.8684 

1.28947 

0.0260843 

3.0192 

500 

95.0356 

0.581249 

0.0196577 

2.9729 

500 

94.7242 

1.07907 

0.0107542 

3.03627 

600 

96.3209 

0.530806 

-0.0227588 

2.92353 

600 

96.0485 

0.957456 

0.0285771 

2.98779 

700 

97.2547 

0.487275 

-0.00846606 

2.98083 

700 

97.0346 

0.85989 

0.0016651 

2.98035 

1000 

99.1397 

0.356686 

0.0636354 

3.13519 

1000 

98.9859 

0.643946 

0.0312824 

2.952 


(c) 


n = 2N 



(d) 


Figure 12: Well-conditioned matrices. Numerical calculations when n = 2 TV. All calculations are 
taken over 32, 000 samples, (a) Histograms for the simulated condition number of LUE. (b) Histograms for 
T e E,N. 2 N for e = 10~ 14 for LUE. The discrete nature of the distribution persists in the N — ► oo limit, (c) 
A table of the computed mean, variance, skewness and kurtosis of T c ^e,n, 2 N for both LUE and PBE. This 
table indicates that the mean is bounded and the variance is monotonically decreasing as a function of N. 
(d) Histograms for t C! e,n, 2 N for both LUE and PBE plotted on the same axes for N = 400. From this plot 
it is clear that due to the discrete nature of the distributions they will not coincide as in Figure [To| 
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A.3 Critically-scaled random matrices 


Finally, we consider the distribution of T e E N jv+i VIlvj > *- e - n = N + [VAcN\ (c = 1) in the case of LUE 
and PBE. In Figure 13(a) we examine the condition number which we know by Theorem 1.3 has Tracy- 
Widom fluctuations for LUE. In Figure 13(b) we examine the distribution for X) E N jv+Lx/Uvi f° r LUE. The 
calculations show a limiting form for the halting time (see also Figure 10). In this case, the moments of 
the condition number are unbounded as N — > oo but the limiting distribution is not heavy-tailed and a 
non-trivial limit exists. Numerical experiments also indicate that this phenomenon persists for the scalings 


n = N + L(4clV) 7 J if 0 < 7 < 1. Universality is also apparent from Figure 13(c) 



LUE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

PBE 

N 

Mean 

Variance 

Skewness 

Kurtosis 

100 

133.573 

1.16093 

0.0163642 

3.09041 

100 

132.235 

2.50567 

-0.276531 

3.40614 

200 

258.119 

16.21 

-0.0752041 

3.02723 

200 

255.512 

30.5358 

-0.0575893 

2.9452 

300 

359.573 

33.8961 

-0.00898283 

2.96971 

300 

356.774 

67.6515 

0.0169258 

2.99919 

400 

442.731 

53.6034 

0.00574321 

3.04426 

400 

440.032 

106.758 

0.00737936 

3.00197 

500 

525.601 

75.9278 

0.0095558 

2.96023 

500 

522.989 

146.83 

0.0309066 

3.00585 

600 

599.448 

95.4292 

0.0288976 

2.99469 

600 

596.929 

192.582 

0.025091 

2.9947 

700 

665.955 

119.276 

0.0245652 

3.01978 

700 

663.378 

230.798 

0.0233565 

3.01858 

1000 

838.013 

175.635 

0.0299194 

2.98757 

1000 

835.867 

348.123 

0.0497676 

2.99017 


(c) 


Figure 13: Critically scaled matrices. Numerical calculations when n = N + \ \/AN \. All calculations 
are taken over at least 32, 000 samples, (a) Histograms for the simulated condition number of LUE. (b) 
Histograms for T e E N Ar+ ^ v 7 jvj for e = 10 ” 14 for LUE. (c) A table of the computed mean, variance, skewness 
and kurtosis of T e E N jv+lvTjvj f° r both LUE and PBE. This table indicates that the mean is unbounded and 
the variance is monotonically increasing as a function of N. Furthermore, the kurtosis is close to 3 for both 
LUE and PBE. This is yet stronger evidence for the universality observed in Figure [lOj The computations 
for N = 1000 use 64, 000 samples. 
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B The Parametrices 

In this appendix we present the asymptotic calculations for the Airy and Bessel parametrices. 


B.l The Airy parametrix 

We use the connection formulas (ui = e 27 "/ 3 ) and refer to Figure [ 3 ] for the sectors 1,11,111 and IV: 

Ai(£) = -wAi(w() — uj 2 Ai(w 2 £), 

Ai'(£) = —w 2 Ai'(u;£) — wAi'(u; 2 £). 


Here Ai(£) is the Airy function. We also calculate the large-£ asymptotics of the matrix functions in (3.1) 
using 


Ai(£) = —= C 1/4 e~** 


2 , 3/2 


2^/tt ' 
1 


Ai , (0 = -^ 1/4 e"* f 


2 , 3/2 


2 y/n' 


(1 + o(r 3/2 ) 

(i + o(r 3/a ) 


(B.l) 

(B.2) 


uniformly for | arg£| < tt — e for any e > 0. 
• For £ £ I we find 


Pm(0 = 


_J_y-i/4 p -§« 3/2 1/4 


2 t3/2 

3 s 


es 


u -nH + 0({- a/2 )) ■ 


Note that w 2 £ £ IV and we must keep with the convention that arg£ £ (—7r,7r]. We used (w 2 £) 3 / 2 = 

|£|3/2 e f (argj —277/3) = _£3/2 &nd ^2^±l/4 = |£| ±l/4g± f (arg £-2tt/3) = £±1/4^1. 


For £ £ IV we find that w£ = |£|e 1 ( ars ^+ 2 s r ) £ I and 

-v ’e --' c3/2 

T Ai (£) = 




2y/n 


1/4 —l/4 e |« 3/2 


2^ 

,„ 1 / 4 . 


- 2 ^ 1 / 4 e -§£ 3 / 2 I ^ C 1 / 4 e ^ 3/2 


3/4 (1 + o{c 3/2 )) • 


This follows because (w£) 3 / 2 = |£|3/2 e i(3/2arg£+,r)) = —£3/2 and ( w £)± 1/4 = |£|3/2 e ±i(l/4arg|+7r/6)) = 
£ ±1 / 4 e ±i ^. 

• For £ £ II we rewrite the matrix using the connection formula because the asymptotics for Ai(£) are 
not uniformly valid here: 


T Ai (£) = 


uj 


-cr 3 /4 


1 0 

-1 1 


Ai(£) Ai(w 2 £) 

Ai'(£) w 2 Ai'(w 2 £) 

— u> Ai(w£) — ui 2 Ai(w 2 £) Ai(w 2 £) 

— ui 2 Ai , (w£) — ui Ai , (w 2 £) w 2 Ai , (o; 2 £) W 

— Ai(w£) — Ai(w 2 £) 

—uj Ai'(w£ ) —uj 2 Ai'(w 2 £) 

— Ai(w£) — Ai(w 2 £) 

—uj Ai'(w£ ) —uj 2 Ai ' (w 2 £) 

—wAi(w£) Ai(w 2 £) 


uj 0 
2 -1 


UJ 


UJ 

1/2 


-0-3/4 

W - CT3/4 

0 


1 0 

-1 1 

1 0 

-1 1 


w 1 / 4 
0 —w -1 / 4 


,-<^ 3/4 


—uj 2 Ai'(w£) uj 2 Ai'(w 2 £) 

For £ £ II we write w£ = |£|e larg £ _4l7r / 3 £ IV and w 2 £ = |£|e larg ^ _2l7r / 3 £ I to compute 


T Ai (£) = 


^J_£-l/4 e -§t 3/2 


2^/tT' 


1/4 -l/4 e |« 3/2 


2 v^7r ■ 

a 1 / 4 . 


^ 2 ^ C 1 / 4 e - i £ 3 / 2 ^£ 4 / 4 e §£ 3 / 2 


3/4 (1 + o(c 3/2 )) • 
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• The calculation for £ € III is similar to this: 


Pm(0 


Ai(0 

Ai'(0 


—co 2 Ai(w£) 
- Ai'(wC) 


CO 


-cr 3 /4 


1 0 
1 1 


—co Ai(w£) — co 2 Ai(w 2 £) 
—to 2 Ai^cj^) — co Ai^w 2 ^) 

— Ai(w£) — Ai(w 2 £) 
-wAi'(w^) —co 2 A\(co 2 ^) 


—co 2 Ai(w£) 

oo- a 3/4 

' 1 

0 ' 


Ai (wf) 

1 

1 


co co 2 

u~ a 3/4 

' 1 

0 ' 


co 2 0 

1 

1 


— Ai(w£) — Ai(w 2 £) 
-wAi'(w^) —co 2 A\(co 2 ^) 

—co 2 Ai(w 2 £) —co 2 Ai(w£) 
—wAi'(w 2 £) — Ai'(w£) 


0 1 
1 0 


CJ 2 C0 


-0-3/4 


,— 0 - 3/4 


For £ £ III we write = |^|e 1 arg €+ 2l7r / 3 g IV and w 2 £ = |£| e lar g£+4i7r/3 g I to compute 


Pm(0 


-^e 1 / 4 e-i ^ 3/2 f ^ 1 ^ 3 / 2 


C0 ~ a3 /4 (i + 0(C 3/2 )) ■ 


B.2 The Bessel parametrix 


We refer to Figure [4] for the sectors in the complex plane for the Bessel parametrix. To state the precise 
Riemann-Hilbert problem (Ricmann Hilbert problem |3.2| ) that is solved by Pnessei in (3.2) we have to address 
asymptotics as a —t 00 . For any e > 0 (see [20]) 


I a(az) 
l' a (az) 
K a (az) 

K ' a (az) 

V(z) 


1/2 


9 arj 


(i) (TTiy7i (1+0( "' 1)) ' i"s4<y-«. 


( 1+-- 2 ) 1 '- 1 
(3 


/ 1 \ 1 / 2 pQt7 1 

I - ' — (l + O^- 1 )), | arg z\ < -7r — e, 


\27ra ) 

1/2 e 


— CiT) O 

(l + O^- 1 )), | arg z\ < -7r — e, 


2 a / (1 + z 2 ) 1 / 4 

1/2 e ~ a v 


-< 1+ ‘ 2 > ,/4 (£) 

(l + z 2 ) 1/2 + log 


(1 + 0(z i )), | arg z\ < -7r-e, 

z 2 


1 + (1 + 2 2 ) 1/2 ' 


(B.3) 


1 /q 

Here the branch cut for (l + z 2 ) is taken on the segment [—i, i]. We compute the asymptotics of each 
factor in the above expansions replacing z with M z/a. for M = iV+|(a+l). Assuming |z| > 5 > 0 with 
I argz| <\i t — e and M/a —> 00 (see {B.31), we find the following: 
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To further simplify things here, use the scaling so that M/a = O(a). From this we have 

2\ V4 


1 + 


1 + 


o\ 1/4 

(Mzy\ M 


-z 1 / 2 (l + O (a -2 )), 


M 


= —z (l + O (a 2 )) 
a 

2c 


arj —2 ) = Mz -b 0(a ). 


Thus, if 1 2 1 > S > 0 we may simplify the expansions 

1/2 Mz 


/ i \ 1/2 Mz 

la(Mz) = Ua(Mz/a)) = J ^~ 2c ' z {l + ©(a" 1 )), | arg 

1 \ 1/2 e Mz 


z\ <2 n ~ e ’ 


l' a (Mz) = l' a {a{Mz/a )) = 


2t tM . 


^,i/ 2 e_2c/z ( 1 + 0(a~ 1 )), | arg 2 1 < -n — e, 


K a (Mz) = K a (a(Mz/a)) = 


= (JL\ 


1/2 e 


-Mz 


V2 M) 


—-e 2c/z (l + C>(a- 1 )), | argz| < 7r — 


/ 7t \ i/ 2 e 


= K(a(Mz/a)) = - (^) 


-Mz 


yl/2 


s 2c /*(l + 0{a~ 1 )), |argz|<7r-e. 


Similar calculations must be completed for the Hankel functions. The simplest way to obtain expansions, is 
to use the fact that for | argz| < 7 t/3 


H«(*) = -e- Q7d/2 K a (ze~ iv/2 ), 

7T1 

H«'(z) = - 2 e- a7Ti/2 K' a {ze- iv/2 ), 

7r 

H V>{z) = ~e a * i l 2 K a {zj*l 2 ), 

7T1 

H f(z) = - 2 e a ^ 2 K ’ a {ze™/ 2 ). 

7r 

Then, | arg(±zi)| < 7r/3 + 7r/2 < ir and we can use the expansions for K Q above as they are still uniformly 
valid. We are ready to compute the asymptotics of Feessei with a both tending to infinity, scaling £ in 
appropriate way. We compute: 

• For £ £ I U IV we have 

i. (^l_) 1/2 £—l/4 e —2M£ 1/2 e c{- 1/2 

2 (nM) 1 ' 2 e /A e- 2M Z 1/2 e*- 1/2 


PBesse\{M 2 £) = 


1 /_J_t - 1 / 2 t-l/4 p 2M$ 

2 UmI e 


i ( 7 rM) i/2 £ 4 / 4 e 


1 / 2 e _cr 1/2 

1 2M{ 1/2 e -c5“ 1/2 


x (/ + 0 (a- 1 )). 


• For £ £ II we have (e 17r / 2 (—^) 1 / 2 ) V 2 = e 17r / 4 (£e 17r ) 4 / 4 = £ i/ 4 e 17r /2 = i/4 anc j 


PBesseliM 2 ^ 

l e -«i/2 ^_L_^ 4 ^ 2 ^—1/4g2M£ 1 / 2 e —c£ -1/2 ^ e «7Ti/2 1/2 ^-l/4g-2M{ 1 / 2 e cr 1/2 

ie -a7Ti/2 (, t M) 1/2 ^l/4 e 2M « 1/2 e - c r 1/2 e ™/2 (ttA f) 1/2 ^l/4 e -2M{ 1 / 2 e cr 1/2 

e |a^ 3 x (J + 0( a -l)). 
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• For £ £ III we have (e 17r / 2 (—^) 1 / 2 ) 1 / 2 = e 17r / 4 (£e 1 ’ 1 ’) 1 / 4 = £ 1 / 4 e 17r / 2 = i£ 4 / 4 and 


PBe SS el(M 2 0 

l e ani/2 (:-l/4 e 2M£ 1/2 e -cr 1/2 

ie^vri/ 2 ( 7r M) 1/2 ^ 1 / 4 e 2M « 1/2 e" crl/2 

x (/ + 0( a -l)). 


| e -a7ri/2 ^j 1 / 2 ^-l/4 e -2M5 1/2 e cr 1/2 

e ~a7ri/2 ( 7rA f) 1 /2 1/2 


B.3 Proof of Lemma 13.II 

We consider 

P Ai (M 2 / 3 f^(z)) = f^(z))~^E M {M 2 / 3 fUz))e- M ^ (z)a3 , 

ZyJ 7T 

and seek a function M A \(z) so that 

S<_(*) := M Ai (z)P Ai (M 2 / 3 f^(z))e^ z)a3 e N ^ {z)rT3 


matches with the outer solution S^. We make the ansatz 
M Ai (z) := 2^D^ 0 1 Af(z)(ip^(z))~ ,T3 


a r 1 / 4 a// 2 


1 “I 


(M 2 / 3 /^(z))i CT3 . 


Using M = N + + 1), (2.17),(2.18) and (2.26) we find uniformly for |z — 1| = S 


£-(*) = D^M(z)(^(z))-° 3 


w- 1 / 4 w 1 / 2 

-W 1 / 4 w 1 / 2 


E Ai {M 2/3 f^(z))e-^ a+1) ^ iz),T3 e^ z)a3 


= D _1 _B Ai (2;)A/’(0)e^ _los ^^^ _ 2( Q + 1 )'! i <-( z )+(“+ 1 ) z -5 Qlo s z ) <T 3 

= D^E^D^S^), 


E Ai (z) = J\f(z){'tp i -(z)) 173 
= I + 


a;" 1 / 4 w 1 / 2 

-W 1 / 4 w 1 / 2 


Em{M 2 ^ 3 {z)Y 3 N~ 1 {z) 


This calculation depends critically on (2.261. 


We now show that M A i(z) is analytic in a neighborhood of 2 = 1. First note that 

• •-V 4 


W 1 / 4 w 1 / 2 

-W 1 / 4 a; 1 / 2 


0 1 

-1 0 


1 -i 
1 i 


Then by direct calculation 

Af(z)(tp^(z))~ a3 
From this we find that 

M m (z) = 2 V W /4 D- 1 


' 1 

—i 

-1 

0 

1 ' 


' -5 

1 

i 



1 

O 


. K z ~ I) i . 


- I'l'l-CTs/ 4 


(z(z- 1 )) 


K z ~h) i 


/CfO, 3/6 /" 

\z(z~ i) 


0 - 3/4 


(B.4) 

(B.5) 


Then, it remains to show that the ratio is analytic and does not vanish in a neighborhood of 0 = 1. 

But this follows directly from the expression for f<~(z) in terms of a power series. 
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Finally, we check the jumps of S < -(z). In the following calculations we leave of ± signs for boundary 
values for functions that are analytic in a neighborhood of the point under consideration. Recall that the 
contours F| and Fj. in a neighborhood of z = 1 are defined in FigurejhJ For z £ T-j-, /<_(z) £ 72 and therefore 

S+(z) = M Ai (z)P£(M 2 / 3 f^(z))e^ (z),73 e N ‘ f ’^ (z) ' T3 


= M Ai (z)P^(M 2 / 3 U(z)) 


1 0 
1 1 


e 2 


iw(z)a 3 Nc/>^(z)cr 3 


= M Ai {z)P^ 1 {M 2 / 3 f^(z))e^’ iz) ' T3 e N ^^ z),T3 e-^ z)rT3 e- N ^^' 73 
= S^L(z) 


1 0 
1 1 


1 0 

e 2 Nifi^(z)+w(z) y 


The same calculation follows for z £ Tj,. For z £ (1,1 + 5) we have /«_(z) > 0 so 
S+(z) = M Ai {z)P£ i {M 2 / 3 f^{z))e^ z > 3 e N< ‘ , ^ ,T3 


= M Ai {z)P^M 2 / 3 fUz)) 


1 1 
0 1 


e 2 


^w(z)<T 3 N 4 ><- (2)0-3 


= M Ai (z)P^{M 2 / 3 f^(z)]e^ z)rT3 e N,l, ^ z) ' J3 e-^ z) ' T3 e- N ‘ l> ^ z)iy3 

= S^{z) 


1 1 
0 1 


I e 

0 


- 2 N<fr<-(z)-w(z) 

1 


Finally, for z £ (1 — 5, 1), /<_(z) < 0 and using that </>+(z) = —(/)<_{z) 
S+{z) = M Ai (z)P^(M 2/3 f^(z))e^ z)a3 e N ^ z),T3 
= M Ai (z)P^(M 2 / 3 U(z)) 

= M Ai (z)PA i (M 2/3 /^(z))e5 rii ( 2) ^e iV ^( z ) CT 3e-|*W^3e^tW^ 

0 


0 1 

-1 0 


e ±w(z)a 3 e Nc/>+_ (2)0-3 


= S<_(z) 


w(z 


e ~w(z) 

0 


e \w(z)cr 3 e N<j>^{z)a 3 


e 2 


±w(z)cr 3 e N 4 ><-(z)<j 3 


0 1 

-1 0 


e \w(z)a 3 & N <t>^(z)cr 3 


We note that both S and are analytic in B{\, S) \ TAi and are continuous up to the contour. This 
follows from the fact that has unit determinant. To see this, note that P Al has constant determinant 
in each sector of C \ Eai by Liouville’s formula and the fact that this constant is the same in each follows 
from the fact the jump matrices have unit determinant. Thus S<-(z) has a determinant that is independent 
of both M and 2 and the determinant is found to be unity by examining its asymptotics. 


B.4 Proof of Lemma 13.21 

We consider 


and 


■Mz) = M Beaael (z)P Besse ,(M 2 f^(z))ei^ 3 e N ^ z ^ 3 , 

M Bessel (z) := D-^izX^iz))-* 3 | l (M 2 U(z))^ 3 i r^ 3 , 

w(z) := (2a + 2)z — a log_^ z + (a + l) 7 ri, 


PBessel(M 2 0 
^Bessel {M 2 £) 


(nM)-i* 3 C ia3 EBe SS ei(M 2 0e 2M S 1/2 ° 3 e- c t- 1/2 ° 3 , 
■i(l + 0 (a- 1 )) i(l + 0(a~i))- 
_ i(l + 0 (a“ 1 )) l + Oia- 1 ) / 
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We first check the asymptotic behavior for |^| bounded away from zero. 


S^(z) = 


1 i 

2 2 

i 1 


E Bessel (M 2 /_>(*)) 


x e ( _ 5( a + 1 )^->( 2 ) + ( a + 1 ) 2_ 5“ 1 °g->( z )+5(“+ 1 )' 7ri ) o '3 e 2c/0_ > (i)(T3 

= D- 1 £JBessei(2)AA(z)(V'^(-)) _<T3 e ( "5 (a+1) ^ (z)+(a+1)i: “5 alos ^ (z)+ 5 (a+1)7ri),T3 e 2c/0 - (z)<T3 

= D- 1 ^Bessel (z)D 00 S 00 ( Z )e 2c ^ 


E Besse[ {z) = Af(z)(l/j^(z)) CT 3 2 2 


EBessel(M 2 (z))(lp^ (z)) a3 JV~ X (z) =7 + 0(a" 1 ). 

This follows from (2.25). Note the extra factor of Q 2c /'t‘^( z ') cr 3 when comparing this with (B.4). 


We now check that M Besse i(z) is analytic. So, we must consider 


N{z){^{z))-° 3 


0-3 

[ 1 -11 


1 0 ‘ 


l-i 1 1 


i(l - 2 z) 1 


(z(z- I))-* 3 / 4 . 


From this, the question of analyticity of M Besse i(z) is reduced to the question of analyticity of the function 


M*) 

z(z-l) 


1/4 


This is clearly analytic in a neighborhood of z = 0 because /_>(0) = 0 but /(_>(0) < 0. 

Now, we check the jumps. Recall that T-j- and T| in a neighborhood of z = 0 are defined in Figure [8j For 
z £ T|, f->(z) £ /?3 (see Figure [4]). The limit to from above (+ side) is the same as limit into fa from 
below (— side) so that 


—e 


1 0 

— OTK\ Y 


S+(Z) = MBes S el(^P B e S sel(^ 2 A(s))e^^)" 3 e JV ^W CT3 
= M Bessel (z)P+ ess JM 2 U(z)) 

= Sl>(z) 


—e 


1 0 

— a- 7 ri 


e ±w(z)a 3 e N (2)0-3 


e ±w(z)cr 3 e Nc/>^(z)cr 3 


1 0 

g— (a+l)iri+ 2 Ntf)->(z)+w(z) 


Then we find that e ( a + 1 )' n ’ 1 + 2Nr f’^( z )+ w ( z ) = e 2 N<f>->(z)+u>(z) and ^g j um p agrees with the corresponding 
jump for S(z). For 2 £ Tj., f~^{z) £ /3i and 

M*) = M Bessel (z)P~ essel (M 2 U (z))e^^ 3 e N ^^ 3 


Bessel v 

= -^Bessel(^)FBe Sse i(M f^(z)) 


= (£) e _ 2 ,& ( z )' T 3 e _Ar '^ z ) <T3 


1 0 


—e 


e ±w(z)cr 3 Q N rj>^(z)a 3 


-e 


1 0 

OtTVl 2 


e 2 


\w{z)a 3 cj>^{z)cr 3 


1 0 

g(a+l)7ri+2AT0_>(2!)+uj(z) ^ 


= M*) 

Then because log^ 2 = log^ 2 + 27ri we have e (“+ 1 )^+2 n^{z)+w{z) = e 2 JV<M*)+<&(*) and agaill) this is t h e 
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same as the corresponding jump for S(z). Finally, for z G (0, S) we have f-*(z) € /?2 and 




= M Bessel (*)P+ (M 2 A( 2 )) 


~lw-(z)a 3 - N 4>_ > (z)a 3 


= S^(z)e 2 —v-/“* e 


0 -1 

1 0 

0 -1 

1 0 


e 2 


i iii + (z)ct 3 p iV<£+ (z)a - 3 


e 2 


i«; + (2)cr3 e Ar0+(2)0-3 


= ^(~) 


0 e ~™( z ) 

-e™( z ) 0 


This follows from <j>%(z) + <j>^(z ) = 0 for 0 < z < 1 and ^(w + (z) + w-(z)) = w(z) + 7ri. It also follows that 
det S-t(z) = 0, see |2Q1 Sections 10.5 and 10.28]. Finally, we need to check that S-^(z) satisfies: 

• S-t(z) = 0(1) as z —> 0 from outside the lens and 


S^(z) 


1 0 

± e 2N(t>^(z)+w(z) j 


= 0(1) as z —> 0 inside the lens. The (+) sign is taken for z in the 
region enclosed by [0,1] and T^ and the (—) sign is taken in the region enclosed by [0,1] and T^. 

From outside the lens, we have for z G P(0, S) \ (0, S) 


S^-(z) = M BesS e\{z) O 


\ z \»/2 |.|-a/2 

\ z \ a /2 U|—a/2 


X O 


|-«/ 2 0 


0 


|a/2 


^ e (( a + 1 ) z + N <t>^(z))a 3 _ Q(\y 


The statement inside the lens follows from the fact that S->(z) 
tinuation of the function defined outside the lens. 


1 0 

± e 2N</>-y(z)+w(z) -y 


is the analytic con- 


C Estimates at the hard edge 

The following facts are of use below. As a —>• 00 we use 


4C 


1/4 


' 1/3 (Ai(a 2 / 3 C) + 0(a" 4/3 )) , t > 0, 


J a (at) = l —^ I a 
3a (at) = ~ (1^2) a_2/3 ( Ai, ( a2/3 0 + 0(a _2/3 )) , t> 0, 




|C 3/2 = / v ^ " da, 0 < t < 1, 
3 Jt s 

4 Vs 2 ^ 7 ! 


"(-C) 3/2 = / —- ~ds, t> 1 

3 Ji s 


(C.l) 


This expansion is uniform for t G (0,00). In preparation for the proofs of Lemmas 4.1 and 4.2 we look to 


obtain global estimates on V and W. Assume z G (0,5') where S' < S < 1/2. The following estimates are 
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straightforward 


—i0i>(z) = y-ds — 4v^, 0 < z < S', 

—i > A\fz\l — <5| 1//2 > 2y/z, 0 < z < S', 

| C 3/2 = f o<t<i, 

3 J t s 3 

C>2 1/3 (l-t), 0 < i < 1, 

^C 3 / 2 > [ ^ S ds = — log2f + log 2 + 2 log(l + Vl — t) — 2-\/l — t, 0 <t< 

3 Jt s 

/Q\ 2 / 3 

C>(j] | log 2f| 2 / 3 , 0 < t < 1. 

It follows from ES Section 9.7] that there exists a constant C > 0 such that 


i(i+M) 1/4 Ai(z)i<c{ 
id + n)-/‘a,' ( x)i<c{ 

Define 

Mt) = (r^) V4 Ai (« 2/3 0- ^(t) = 1/4 Ai'(a 2 / 3 C). 


Let C > 0 and t < C. For t £ [1 — e, C], Aj and A 2 are bounded as follows 


Ai(t) < C max_ 

tG [1 — £,0] 


4C 

1 -t 2 


1/4 


A 2 (f) < C(1 + |a 2 ^ 3 C|) 1//4 max_ 

£G [ 1—e, 


4C 

1 -t 2 



For t € (0,1 — e) and a constant C e > 0, 

|Ai(t)| < y/2(l - t 2 )- 1/4 |C| 1/4 | Ai(a 2 / 3 C)| < C e e~ a ^ 3/2 < C c e~ a ^ 3/2 

iA 2 (t)i<^( i -t 2 ) i/ 4 ici- i/ 4 i Ai ' < on)i 

< 71 (1 - t2)1,4 (l + |aC»;0-i/l (1 + l« 2,3 fl) _1/4 l Ai '(“ 2/3 OI < &e-K 


1^3/2 


because 


ICI 


-1/4 


< ce a z 


c > 0. 


(1 + |a 2 / 3 C|) _1/4 

In summary, for any > 0, there exists a constant C = C'(('i) > 0 such that 

\Ai(t)\ <c{ e _a 3 |c|3/2 , if C > 0, 

l ' 4l(t)l - l i, if -(,«-w < C < o, 

\ a (f\\ < c I e_ “® IC|3/2 f if C > 0, 

|A2(i)| - C i i, if -Ci«- 2/3 <C<o. 

Note that if £ > — £ia _2//3 then t is bounded, i.e. t < C. The constant £i is fixed below (C.9). 
C > 0 then |Aj(i)| < Ce _ “3^l 3 2 < C so that the second inequalities above hold for — £ia -2 / 3 
remarks apply in (C.5| below. 


(C.2) 

1, 


Note that if 
< £. Similar 
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Define f as a function of 2 by t(z) = —yfi <!>%.(z) so that 3 a (at) = Jq,(— iM(/>+ (z)) to match (|4.3[) and 


(4.4). From the estimates in (C.2) it follows that t(z) < 4so that if t(z*) = 1 then (jfj) < 2* and 
( \{t(z )) > 0 for 2 < z* , as t(z) is an increasing function of z. Furthermore, we know that £ > 2 1 / 3 (1 — t ) for 
0 < t < 1 so that 


COO)) > 2 1/3 ^1 - , for 2 < (£j^) and 


-Cic*- 2/3 < COO)) if - Ci a_2/3 < 1 - —vO. 


-2/3 


4 M 


For 4 —Jz < 1 we obtain 

Oi V — 


3 -«| ici 3/2 < e -« K 1 - 4 “^) 7 . 


From the last line of (C.2) we have 

e —“ 3 IC | 3/2 < e -ia|log 2 t| < ^ t) a / 2 < e"^ 108 ^^)^]. 


(C.3) 


(C.4) 


As we will see (C.4) is useful near z = 0 and ( |C.3 ) is useful for slightly larger values of z. Using ( C.l| ) we 
have for a constant C\ and any 0 < C 2 < (yyf) 

I Ja(-iM<^+ 0))l = a -1 / 3 |A 1 (t)||l + 0(a- 4 / 3 )| 

e“5 1 °g[(“)v^] i if 0 < z < C 2 , 


<Cra V 3 <| e-^0- 4 ^^) 372 , if c 2 < z < (^) 2 , 

1, if -Cia~ 2 / 3 <l-^-^z<0, 

0+0) J' a (-iM0+0))| = 2M- 1 a 1 / 3 |A 2 (t)||l + 0{a~^)\ 

e “Ooif 0 < z < C 2 , 

<CrM-V/3<( e -aI(l-4M^)3/^ ZCzKzK^) 2 , 

L, if — CiCt -2 / 3 < 1 — ^zr\[z < 0, 


(C.5) 


and C 2 is determined below in (C.10). 


In what follows, we also need an estimate on W', see (4.4). So, we consider 
W\z) = 


-iM0+' 0)J0 (-iM<£+(z)) 


7T0+ 7 (z)3a (—iM (f)%{z)) + CO/OX - iirM<f>t>(z))J a "(— iM0+ (z)) 
It also follows that zJ a "(z) + J a '{z) = z~ x (a 2 — z 2 ) J a (z) and therefore 


e (iV+|)7ri^ 1 /2_ 


W’(z) = 


hr M 


~ iM l+ ri* < A0O)Ja / (-iM^>+ (z)) 

<piO) 

A,f2/4,+ /'^\2 


^iO) 


(a 2 + M 2 {<j)\{z)Y) J a (—iA/O+O)) 


eC-^+D^Tl^ 1 / 2 . 


Then we have for 0 < 2 < 5' 


<o0o) 


</>00) 


<2; 1 \/i — z < z 1 , 
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so that 




< C.a 1 ' 3 


3 ailog[(SM)v^]-l°gz if 0 < 2 < C 2 


e z ^l\ a /» j - 5 II U ^ ^ 02 

x { e-^K 1 - 4 "^) 372 -!^^ ifC 2 <z <(^' 2 


2 =i * ^ v4m; > 
if - Cia- 2/3 < 1 - ^Vz < 0, 


M~ x ^~! [ Z } (a 2 + M 2 (0+(z)) 2 ) J a (-iAf^+(«)) 




< C\M~ x aT x ^\o? + M 2 (0+(z)) 2 | 

if 0 < z < C 2 , 


e -«i(i- 48 fVi) --iogz ; if C2 < 2 < (jfij) 2 , 

if -Cia" 2 / 3 < 1- 0. 


v-1 


We write B(y) = B(x) + (x — y) B'(tx + (1 — t)y)dt and then 

/C N (x,y) = _J_ v ( x ) w (y) _ -}-V{x) ( [ B~ x {x)B'(tx + (1 - %)dt) W(y). 


2ni x — y 2m 


Recall the scaling operator 


c 2 c 2 2 2 / 3 c 




x = 


= — l+x - 


a* or 



and we make the transformation 


)Cn(x, y)dy -> IC N (x , y)dy = IC N (x , y)dy. 


We find 


1 c 2 /2' 2/3 


W«)di = -i«d 9 - . 

Z7ri x — y z7Ti a z \ a 

C.l Proof of Lemma 14.21 


V{x)( / B~ x {x)B'(tx + {l-t)y)dt )W(y)dy. (C.6) 


We prove Lemma 4.2 first. For x in a compact set Q, as IV — > 00 using = 1 + 0(a 3 ) we have 


2/3 


-i<£+(x) =4-\ l + x ( - (1 + O(o~ 2 )) =4-(l+a;a" 2/3 2- 1/3 + e)(a" 4/3 )), 


—i— 4>%{x) = (1 + 0(a 1 ))(l + a;a 2 ^ 3 2 1//3 + 0(a 4 ^ 3 ) = 1 + xa 2 ^ 3 2 4 / 3 + 0(a 4 ). 
a 

Then, £ = 2 1 / 3 (1 — i)(l + 0(i — 1)) so that 

£(i( x)) = —a~ 2 / 3 x + 0(a~ x ), 

f 4 C0(a)) \ 1/4 


\ 1 — t 2 {x) J 


= 2 1/3 (1 + C>(a- 2/3 )). 


Thus by (C.l) 


/ 9 \ 

J a (—iM^ (£’)) = J a (at(x)) = ( — J Ai(—x) + ©(aC 2 / 3 ), 

/ o \ 2 / 3 

7T(/>+ (x)J a , (-iM^+ (£■)) = 7ri-^:t(x) J a ^(crt(x)) = -7ri^r ( - j (Ai'(-x) + £>(a _1/3 )). 
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Here convergence is uniform in x for a; in a compact set. The following estimate also follows for — (^)~^ 3 < 
x < L (recall (3.51) 


a 2 + M\4>%{x)f = a 2 - 




2 \ 


2/3> 


(l + 0(a- 2 )) 2 . 


Then from = 1 + ©(a" 1 ) we have 


a 2 +M 2 (</>+(x)) 2 =a* -a\l + 0{a~ 1 )) ( 1 + s ( - 


2/3> 


a 2 + M 2 (</>+ (x)) 2 = a 2 - a 2 - xa 4/3 2 2 / 3 + 0(a 2 / 3 ) = -a A/3 2 2 / 3 (x + 0(a~ 1/3 )), 
|a 2 + M 2 {(f>Y{x)) 2 \ < Ca 4 / 3 (\x\ + l) 2 , 
for some C > 0. We proceed to calculate 


(C.7) 


^'(x) _ 1 [l f <^ (1 + xa _ 2/3 + 0(a _ 1)} = 1^(1 + ©(a" 2 / 3 )), 


</>+(£) 2 
r <()+I 
r _l£±/(*) 


x c 


2/3 


-iM^|| 0 + ( 4)J Q '(-iM^+ (£)) = - ~ (± ) (Ai'(-x) + G(a- 1 / 3 )) 


It rM 1 ,'7/ (a 2 + M 2 (</>+ (£)) 2 ) J a (-iM^+ (£)) 

•few 

1/3 


= —17T- 


a 

cAM' 


4 / 3 ^ Ai(-*)(* + 0{a~ 1/3 )). 


Combining this discussion we have proved Lemma |4. 2 


C.2 Proof of Lemma 14.11 


Next, we turn to global estimates after the change of variables x. We must consider expression in (C.51 after 
this change. To examine the components of V and W we use 


^ 4 Me < 41V + 2v / 4o N + 4 < 1 + 2(4c7V) _1/2 + 4(4clV)- 1 ^ o 
“ (V4dV - 1)2~ ~ 1 - 2(4ciV)- 1 /2 + (4ciV) _1 < ’ 


v2 - 


(C.8) 


for sufficiently large N to see that 

log 


/ 8M 

V « 


<-!„g 


16 [ 1 + x ( — 
a 


2/3> 


Now, assume — ( a/2 ) 2 / 3 < x < —1 and then for sufficiently large N 

2/3\ 1 / 2 

;) 


M n ; 

1 — 4— Yx = 1 — 


4 Me 


a 


a “ 


1 + a; ( — 
a 


= i-d + d4f /3 ' 


1/2 


a / 


1 - 


^)h(!f) 


> 


(9 


2/3 


1 (a\ 2 / 3 ( lMc\ 


~2 X + 


*<§r 


(I) 

—(iro-^ 5 ) 


■ 




2/3' 


1/2 


1 / 2 ' 
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This follows because 1 — (1 + y ) 1 / 2 > —y /2 for y < 0. Then 1 — AMc/a 2 = 0(iV - 1 / 2 ) and the second term 
in brackets vanishes as fV —>■ oo. Thus, for sufficiently large N 


1 — A—VS; > — I — 


2/3 


X + 1 


> 0. 


Then the inequality 0 < x < (jf/) 2 holds and 


i -M A M 

exp | -a- [ 1 — 4— 


Vi) 3/2 )<e*p(4l* + l| 3/2 )- 


Next, for — 1 < x < L consider the inequality 


9/ , 4 M r- 

-Cia < 1 - Vi. 

a 


The right-hand side is a decreasing function of x so we consider 

2/3 \ 1 / 2 '' 


-Cl < « 2/3 1 - 


4 Me 


1 


( / 2/3 \ 1 /^ 

1 — (1 + 0(7V -1 / 2 )) f 1 + L ) 

( 1+ © 2/3i )' /2 ) +0( °" /S> 


= a 2 / 3 I 1 - 


= -2 - 1/3 L + e>( a - 1/3 ). 


Given L, we choose Ci so that this condition holds and the estimate for — Cic* 1 < C < 0 in (C.5 
used for — 1 < x < L. Then 

f e flog [ 16 ( 1+a; (“) 2/3 )], if - (l) 2 / 3 <x< -&(f) 2/3 , 
\W 1 (x)\ = \V 2 (x)\<C 1 M 1/2 a- 1/3 l e -§|x+i| 3/2 ; if _ 6(|) 2/3 < x < -1, 

[ 1, if — 1 < x < L, 

( e^ log [ 16 ( 1+x (“) 2/3 )], if -(f) 2/3 <z<-6(f) 2/3 , 

|W 2 (x)| = |Vi(x)| < C , i7rM' 1/2 a 1/3 ^ e -||x+i| 3 / 2 ) if - 6 (f ) 2/3 < x < -1, 

[ 1, if — 1 < x < L, 

for any b such that 15/16 < b < 1 by choosing 

C 2 = c 2 cG 2 (1 — b). 


(C.9) 
can be 


(C.10) 


52 





We note that by (C.8) this choice of C 3 satisfies C 2 < (jfy) 2 - We estimate the derivative using (C.7| 


M 


\W[{x)\<^\W 2 {x)\ 


<C 3 M^ 2 a 1/3 % { _L e -|l*+i| 3/ 

c l 1-6 


flog[l6(l+x(f) 2/3 )]-log[l+x(f) 2/3 ] ; if _ (|) 2 /3 < x < _ b ) 


\2/3 


if — 6 (f ) 2/3 < x < -1, 
if — 1 < x < L, 


| W , (i) | s ^ + M^(fc(x))^| |tyi( , ) | 


3 f log[l6(l+x(f) 2/3 )] -log[l+x( J) 2/3 ] 


J, if -(|) 2/3 <x<-6(f) 2/3 , 

<c 3 (M + i) 2 M-v^^ _l_ e -ii,+n 3 / 2 ; if -6(f) 2/3 <x<-i, 

1, if — 1 < x < L, 

for some C 3 > 0. This proves Lemma pl~T| after choosing 6 sufficiently close to unity so that 

(l + | a; |)2 e f lo s[ 16 ( 1+ "(-) 2/3 )]- los [ 1+x (-) 2/3 ] < 1, 

Q /Q Q /Q 

for all a > 4 and — (f)” <£<—&( f) because then W[ and W' 2 have bounds that are independent of 

x. We use 6 = 1 in the estimates for V and W. 

C.3 Proof of Proposition |4.1| 


From (C.6| we have 


iC N (x,y) = -^ V(x)W{y) (-) V{x)(j 1 B-\x)B'{tx + {l-t)y)At\w{y). 

Z7ri x — y Z7Ti or \a J \J 0 J 

We first consider the second term in the expression. We note that /j) S _1 (x)S'(tx+(1 — t)y)dt is a uniformly 
bounded function. Define 


^(*.y) : =2^i«2 («) V ^(J 0 B 1 ( i ) s/ ^ + ( 1_t )^) d ^ w (y)- 


We use 


J e s P+ 1 ! if — 00 < x < —1, 
g[x) = < , ’ — 

1, otherwise. 


By Lemma 4.1 for (x,y) £ (— 00 , L] 2 there is a constant > 0 such that 

\ 2 /3 


1 M f 9\ 2 

|'H(x,y)| < D L a~ 2/3 — ^ ( - ) g{x)g{y) < D L g{x)g{y). 
2ir or \ol) 


For the first term, we assume \x — y\ > 1 and we have 

1 V(x)W(y) 


2m x — y 


< C\g(x)g(y). 


For \x — y\ < 1 we use 


1 V(x)W(y) 
2ni x — y 


J_c 2 _ 
27ri a 2 



V{x)W' (tx + (1 — t)y)dt. 


(C.ll) 


53 
















Next, define t x _ y = tx + (1 — t)y so that t x - y = tx + (1 — t)y. It follows from Lemma 


4.1 


that 


\V 2 (x)W'(t x _ y )\ < C 2 L a 8 / 3 g(x). 


(C.12) 


Next, we consider 


mx)W[{t x -y)\ < C 2 L a 8 ' 3 g(x). 


From this, it follows that there is a new constant Cl > 0 such that 




< C L g{x). 


(C.13) 


Therefore, we compute 


p pL px~ I~1 pL 

/ g(x)dxdy < / g(x)dxdy < 2 / g(x)dx. 

J \x—y\< 1 J — oo J x — 1 J — oo 


The first part of the proposition follows by choosing 


G(x,y) = D L g(x)g(y) + 


C L g(x ), if I* — 2/| < 1, 

G 2 L g(x)g(y ), otherwise. 


Next, for (a;, y) in a compact subset of (—oo, L\ we have from (4.6 1 and (4.7) that V(x) = 0(Ma 2 ^ 3 ) = IT(y) 
so that 


U(x,y) = 0(a 4 ^ 3 ), uniformly in x and y. 


Define t x - y = tx + (1 — f)y and 


£»(«) =-^W™ + o(o- 4 D 

27ti x — y 


where the error term is uniform in x and y. Using (C.ll) we consider We also compute 
c 2 / 2\ 2 / 3 

^2 f - J V(x)W'{t x - y ) = 27ri [Ai'(-cc) Ai ’(-t x - y ) + t x - y Ai(-x) Ai (-t x - v )\ + 0(aT 1/3 ), 

uniformly in x and y. Therefore 

ICa f (x, y) -)• / [Ai'(-x) Ai '(-t x - v ) + t x - y k\(-x) Ai (~t x _ y )\ d t 
Jo 

= Ai(-x) Ai'(-y) - Ai'(-a:) Ai(-y) 
x-y 

uniformly for x and y in a compact set. This proves the proposition. 


D Estimates at the soft edge 


In this section we prove the results presented in Section 4.2 We write B{y) = B(x) + (x — y) f 0 B (tx + (1 — 
t)y)dt and then 


K. N (x,y) = 


V(x)W(y) — 


x-y 




+ V(x)[ B (x)B (tx + (1 - t)y)dt )W(y), (x,y) e (1 - S ,1 + S] 
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Recall that x = 1 + x/(2 2 / 3 M 2 / 3 ). We consider K^{x,y) and we find 

£n(x, y ) = y( ^ (y) + m2/ 3 2 2 /3 ^(*) (£ W(y), (£, y) € (1 - 5,1 + <5) 2 , (D.l) 

for t x - y = tx + (1 — t)y. An alternate expression is 

£*(*,y) = (j[V(Wd*) + jp^nx) ^w\x)B\i x - v ) d*) to, (d.2) 

because V"(x)W(x) = 0. Similar expressions follow with the over line replaced with a tilde: 

JC N {x,y) = M 2P 2 2/3 ^(x) g~ 1 (^)B'(t x - y )d^ W(y), (x,y) € (l + <5,oo) 2 . (D.3) 

Additionally, 


&n(x, y) 


V(x)B \x)B(y)W{y) 

x-y 


(D.4) 


D.l Proof of Lemma 14.31 

For x £ (1 — 5, l + <5) it follows that 


vTT? 


2/3 


(z - 1) < /<-(x) < 




2/3 


0 - 1 ), 


so then for x £ ( L, 5M 2 / 3 2 2 / 3 ), L < 8M 2 ! 3 2 2 / 3 and sufficiently large iV we have 

/ i \ 2 / 3 / i \ 2 / 3 

(vT+l) 1 s " 2,3/ - (i) s (Tra) x - 

We use the following estimates on Airy functions. There exists a constant C > 0 such that 


1, if x < 0, 

_2 3 / 2 

e 3 X , if x > 0, 


I Ai(x)| < C 
|Ai'(x)|<C(l + |x|) 1 / 4 { 


2 : 3/2 


if x < 0. 
if x > 0, 


Then for cj 3 



and a new constant C > 0 we have 


Ai'(M^(*))|<C'(l + M) ,/4 { 

|//(x)Ai'(Af 2 «/^(x))|<C(l + | a ,|) 5 ' 4 { 
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From this we have 


F 2(i )i = F |Fl(i) |< c(1 + wrM -./s{ 

iw;wi<c(i + wrM»/»{^„ mi 
\W 2 m = \fUi)M 2/3 fm*)^(M 2/3 f^m 


< sup 

\ tcG(l—^,l+<5) 


I fU*)\) CcTm^{ X ’ 


xe 


2 3/2 

- —ns v ' 


if x < 0, 
if x > 0, 


Note that 


g-§ c 6 x 3 / 2 Q x _ e -x(lc s x 1/2 -l) 


is bounded and decays as x —> oo and the first part of the lemma follows from these estimates. 
Next we must consider 


W{z) = 
V{z) = 


e h(z)+Ng + (z)+^(N- \vz+ \a log z 

0 

0 e h(z)+Ng + (z)+^t N - ji/z+jalogz 


and obtain the estimate 


h{z) + Ng + (z) + ^jv - m z + log z = -2Nz + !V(21og2 + 1) - (z(z - l)) 1 / 2 - ^a0<-(z) 

< -27V(z - ^(1 + log2)) < —2N(z — 1). 

This proves the last part of the lemma. 


D.2 Proof of Lemma 14.41 


The asymptotics for W and V follow from (4.9). For W we consider 


W'(z) = M 2 / 3 Mb CT3 I" ,,2/3 zf ^T/mV31 ( 

_ M 2 ' 3 f^(z)m(z) Ai(M 2 ' 3 m(z)) 

Then (4.9) along with the fact that M 2 / 3 /{_(±) = 2 2 / 3 M 2 / 3 + 0(M ~ 2 / 3 ) gives the asymptotics for w\x). 


D.3 Proof of Proposition 4.2 


The proof of this follows Proposition 4.1 We first find the function G. Let i?i = [L, dM 2 / 3 2 2 / 3 ] and 
R 2 = ((S/2)M 2 / 3 2 2 / 3 ,oo). Note that R 2 obtained by reducing S. We use 


g(x) := (1 + |z|) 1/4 


if x < 0, 
if x > 0. 


The first case we consider is \x — y\ > 1. We use (D.l) and (D.3) along with Lemma 4.3 and the fact that 
B, B ’, B and B’ are uniformly bounded functions to state that there exists a constant > 0 such that for 
sufficiently large N 


\^n(x,v)\ < D L g(x)g(y), (x,y) £ (i?i x R ± ) U (R 2 x R 2 ), 


(D.5) 
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with \x — y\ > 1. It now suffices to consider (x,y) £ Ri x R 2 , \x — y\ > 1 because an estimate for 
(x,y) £ R .2 x i?i, |a; — y| > 1 follows from symmetry. But estimate (D.5) in this region follows immediately 
from ( D.3| ). It is clear that g £ L l (\L, oo)). 

Now, for \x — y\ < 1 we note that {|x — y\ < 1} C (i?i x I?!) U (i? 2 x 1? 2 ) for sufficiently large AT. Then 
we use (D.2) and (D.3) and find that there exists a constant Gl > 0 such that 

\V(x)w\y)\ < G L M 2/3 {g(x) + g(x)g(y)), \x - y\ < 1, 

and by possibly increasing Dl, we have 

\^n(x,v)\ < D L (g(x)x\ x -v\<i(x,y) + g(x)g(y)) £ T 1 ([T,oo) 2 ), 

where \u is the characteristic function of the set U. Thus G{x,y) = Di,{g{x)x\x-y\<\{x,y) + g{x)g(y)). 

To determine the uniform limit of K-n we use (D.2). From Lemma |4.4| we have 

lim JC N (x,y) = - [ (t x _ y Ai(t x - y ) Ai(x) - Ai'(x) Ai'(t x ^ v ))dt 

N—toc J Q 

where the limit is uniform for (a;, y) in a compact set. Then we see 


- / (■ t x - y Ai (t x - y ) Ai(x) - Ai'(x) Ai'(t x - y ))dt = - 


Ai(x)[Ai'(x) - Ai'( 2 /)] Ai'(x)[Ai(x) - Ai(j/)] 


x-y 

Ai(x) Ai'(y) - Ai(y) Ai'(x) 
x-y 


x-y 


This proves the proposition. 
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